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FOREWORD 


The effort described herein was performed under Contract NAS9-7566. 

The major objectives of this research program were: (l) to experimentally 

identify in "real time" the chemical species present in the combustion 
chamber of a nitrogen tetroxide/aerozine-50 rocket engine during the pre- 
ignition, ignition and post combustion operating periods, and (2) to 
determine, if possible, the cause of severe ignition overpressures sporad- 
ically observed during the vacuum startup of the Apollo RCS engine. 

To achieve these objectives, observation will be' made of the emission 
and absorption spectra associated with specific chemical species present in 
the rocket chamber during the several phases of operation. Specific emphasis 
is to be placed on outlining the physical and chemical processes occurring 
during the preignition and ignition spike periods of operation, and on to 
the extent to which residual propellants or condensed intermediate reaction 
products remain after the engine has been operated in a pulse mode duty 
cycle. 
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I. 


INTRODUCTION 


1.1 Background 

During the development tests of the Apollo vehicle, several hundred 
thousand test firings were performed using The Marquardt Corp. RU-D nitrogen 
tetroxide/Aerozine-50 reaction control system rocket engine. One of the 
problems that occurred during these tests was a sporadic detonation in the 
thrust chamber. This detonation occurred only during high altitude simula- 
tion tests and was not observed during static tests performed at sea level 
firing conditions. In several well documented occurrences these detonations 
resulted in extensive structural damage. 

Even though the basic mechanism responsible for these detonations 
was not fully understood, a number of major modifications of the engine 
were introduced in order to reduce the probability of their occurrence. 

These modifications were based on the hypothesis that unburned propellant 
residues were responsible for the observed phenomena. Such residues were 
observed to form in appreciable quantities during pulse-mode operation 
under simulated high-altitude-start conditions. 

Support of this hypothesis was found in two simple experimental 
observations. The first was the observation that strong ignition spikes 
had never occurred on a clean engine test. This was interpreted to mean 
that the strong ignition spike phenomenon was in some way a result of pre- 
vious engine tests. The second observation was that the frequency of 
detonation was dependent on engine attitude; the frequency increased if the 
engine was tested in a firing attitude such that the propellant residues 
remained in the engine. 
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In the final engine configuration, two design changes were in- 
corporated to reduce the problem. The first was the incorporation of an 
ignition precup in the central portion of the injector plate. The ignition 
precup consisted of an open ended annulus projecting from the axial center 
of the injector plate. - The propellant manifold was designed in such a 
manner that the initial propellant flow started through this annulus. The 
rationale for this design change was that the flow channelization caused 
by this device increased the local stream pressure in the center of the 
thrust chamber during the preignition time period. Because of the expo- 
nential dependence of burning rate on pressure, it was believed that the 
increased pressure would reduce the ignition delay time and promote igni- 
tion near the injector plate, -thus reducing the quantity of propellant 
residues . 

The principal difficulty with this approach was that high-speed 
motion pictures of the ignition process, taken after the production engine 
had been developed, indicated that ignition could have occurred either in 
the precup region or at the nozzle throat. For a throat ignition, it was 
probable that shock heating was involved and that the mechanism of ignition 
was different from that on which precup arguments were based. 

The second major engineering change was the incorporation of ad- 
ditional electrical heating capacity on the injector of the engine. The 
heaters provided sufficient heat to maintain the injector at an elevated 
temperature, apparently removing propellant residuals and engine residues 
by direct volatilization. The test data available to date indicate that 
removal of the engine shut-down propellant residual materials eliminates 
the catastrophic overpressures. However, the use of heaters requires 150 
watts of electrical power for each attitude control engine. Since power 


2 



is at a premium in a space vehicle, this represents a substantial power 
drain on the secondary power source units. 

1.2 Summary of Recent Research 

Although satisfactory engineering solutions were found, the cause 
of the detonation and the fundamental chemical and physical processes 
responsible for the phenomenon remained unknown. As a result, a series 
of research programs were conducted to provide detailed knowledge of the 
chemical and physical processes responsible. 

A n umb er of investigators have provided information about the 
chemical character of engine residues. Seamans and coworkers''' have shown 
that the residues of a MMH/NpO^ propellant combination have a chemical 
composition similar to that outlined in Table I. This elemental composi- 
tion is interpreted to mean that the overall mixture consists of monomethyl 
hydrazine admixed with water and an uncertain nitrate adduct, probably HNO3. 

p 

Using standard Trazuzl block tests, Seamans found that the brisance 
of the residual material was extremely high, and that the explosive force 
was sufficient to completely destroy the test fixtures in some cases. This 
was in marked contrast to results obtained with a well behaved explosive such 
as nitroglycerine. The work of Perlee and coworkers at the Bureau of Mines J 
showed rather conclusively that the infrared spectra of the residue cor- 
responded to that of a hydrazinium nitrate admixed with pure unreacted fuel 

k 

components. In addition, Perlee studied the required shock overpressures 
for the pure samples of several hydrazinium nitrates and found that a min- 
imum shock initiation pressure of 22 kiicbars were required to initiate 
the explosion. 
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TABLE I 


MMH/nTO Residues, o/F = 2.0 


A. Elemental Analysis 


C 16.3 $ 
H 8 . 6 $ 
N 42.8$ 
0 (By difference) 32.3$ 


Empirical Formula 


^ 4 ^ 24 ^ 9 % 


B. Functional Analysis 


Total N-N (as MMH) 

Total ammonic (as CH 2 NH 2 ) 
Total RO 3 (as HNO 3 ) 

Water 

co 2 

wo 2 

53-81$ 

0.44$ 

16.08$ 

15.3$ 

Trace 

Trace 

Total Analysis 

85.6 $ 

C. Derived Composition 

of Residue 

ch 3 nh 3 wo 3 

1.33$ 

MMH • HN0 3 

26.28$ 

MMH • H 2 0 

54.8$ 

MMH 

3 . 22 $ 

$ Total 

85.63$ 
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There is no evidence from the present work that pressures of this 

5 

magnitude were present during the starting transients. Perlee has sug- 
gested that localized pressures of this magnitude may possibly be reached 
via multiple reflections in the chamber. 

Infrared studies performed by Dauerman and coworkers D have suggested 
that, in addition to the hydrazinium nitrate, there is a possibility that 

i 

an azide might be present as a major component in the residues. The 
evidence for this possibility is the observance of a transient peak at a 
wavelength of 4.8 microns, observed under relatively carefully controlled 
laboratory conditions. No evidence for the azide peaks in the infrared 
region was found in the residues examined by either Perlee or Seamans during 
their early work. However, later work by Perlee^ indicated that tetrazine 
derivatives may be present. 

A physical characterization of the residues has been reported by the 
same investigators. The residue material consists of a viscous, yellow 
liquid. Both under high altitude simulation tests and under ambient pres- 
sure conditions, the material is observed to outgas for long periods of 
time. The most probable cause of the gas formation has been assumed to 
involve slow liquid phase reactions between unreacted hydrazine and the 
reaction products. 

Recent work of Saad (private communication) has demonstrated in- 
creased complexity of the residue mixtures. Chemical analysis of the 
mixtures under pseudo combustion conditions indicated the presence of 
dimenthylnitros amine, monomethyl nitrosamine and the probable transitory 
existance of nitrosamine for UTMH, MMH, and N 2 H^ reactions with nitrogen 
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tetroxide respectively. Evidence for these compounds was obtained using 
a combination of infrared, mass spectrometric , and gas chromatographic 
analyses. More stable decomposition products of the unstable intermediates 
were also observed. 

Further rocket testing programs have demonstrated that at least two 
fundamentally different problems are involved in the destructive ignition 
spiking of pulse-mode rocket engines. The problems are associated with the 
processes occurring in the preignition and postcombustion phases of opera- 
tion. Each type of ignition "spike" probably has its origin in a separate 
physical phenomenon. Engine failures attributed to preignition phenomena 
have been observed on "first start" tests performed at the Bureau of Mines 

with a two-dimensional research rocket engine. Because there is no material 

! 

left from previous runs, a failure on a first start must be explained in 
•terms of the combustion preignition kinetics probably involving both gas 
and liquid phase reaction mechanisms. Attention has therefore been focused 
on the gas phase reaction sequences, in part because of the high probability 
that branched chain reaction sequences in the gas phase could lead to the 
initiation of detonation via the formation of unstable chemical intermediates. 
Reactions involving liquid phase components are expected to be relatively 
slow since they are governed by liquid diffusion velocities and/or mixing 
of the fuel and oxidizer components and can play no role in the development 
of high overpressures during the ignition of clean engines. 

Engine failures have been observed which can be attributed to a con- 
densed phase detonation of the liquid residues formed during and after 
engine shutdown. In this case the detonation occurs during the next start up 
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of the pulse mode engine. Such an event may occur by a pressure or tempera- 
ture induced initiation of the "primary explosive" which consists of a mix- 
ture of hydrazine hydrate and hydrazine nitrates and/or nitrites. It- is 
noted that such mixtures are relatively insensitive explosives. In addition 
to a direct initiation of the "primary explosive" residues, a second mechan- 
ism for explaining the liquid phase detonation may be found by assuming the 
\ 

presence of highly sensitive compounds such as hydrazoic acid, trizene and 

hyponitrous acid in sufficient quantity in the liquid phase to initiate 

explosion of the base residues. The present research program indicates that 

(HNO) , (KN) and HN are present in the volatile gases emitted from the 
x x 3 

liquid residues. A recent report by the Bureau of Mines indicates the 
possibility of dimethyl tetrazene as a trace component. If such compounds 

j 

are present as an entrapped gas phase, the interesting possibility exists 
that the liquid phase detonation may involve- the development of a localized 
hot spot during an adiabatic compression caused by engine pressurization. 

There are at least three ways in which hot spots can be developed. 

One involves the adiabatic compression of a gas bubble in a liquid explosive; 

\ : ' j 

a second involves the heat generated by viscous flow; and a third involves 
a mechanical abrasion between crystals. Because only a limited acceleration 
of the fluid can occur during the time period associated with the ignition 
spike, it is probable that this is not the initiation mechanism in this case. 
Since in most of the recent work the residues are present only in a liquid 
form, mechanical abrasions between crystals is ruled out. Because the 
residues are characterized as frothy liquids the adiabatic compression 
mechanism is operative and therefore a detailed analysis of this initiation 
mechanism has been performed. 
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In this analysis it is demonstrated that liquid phase reactions 
produce’ metastable gaseous compounds such as hydrogen azide and triazene 
and create chemically metastable bubbles. The thermal energy increase 
caused by the compression of these bubbles then produces localized "hot 
spots" throughout the liquid. If the temperature of these localized hot 
spots exceeds the minimum activation energy requirements for chemical de- 
composition, then the hot spot undergoes a highly exothermic chemical 
reaction which results in a high energy flux into the relatively inert 
propellant residues. This energy flux can initiate the residues and cause 
catastrophic detonations. Because initiation can be started at several 
locations simultaneously, multiple detonation waves are probably developed 
which can result in a shaped charge effect which would account for the 
characteristically high brisance found in the propellant residues. The 
adiabatic hot spot mechanism is capable of explaining all the phenomena 
observed both under carefully controlled laboratory experiments and during 
hot engine tests. 
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II. 


THERMODYNAMIC ANALYSIS 


2.1 Thermodynamics of Engine Performance 

Propellant oxidant/fuel (0 /f) mixture ratios were optimized with 
respect to selected propulsion parameters such as chamber pressure, area 
ratio, and specific impulse. The theoretically calculated values for 
these system parameters were correlated with actual engine performance 
factors through an empirical correction which accounts for propulsion 
inefficiences and experience on previous engine development programs. 

Two sets of calculations were made using the NASA Thermodynamic 
Analysis Program. The first type of analysis assumed shifting equilibrium 
during an isentropic expansion from engine thrust chamber conditions through 
a DeLaval nozzle. A shifting equilibrium analysis provides a theoretical 
limit to the maximum performance of the system. The two types of analysis, 
shifting equilibrium and frozen equilibrium, comprise a design envelope 
that can be used to estimate the performance and chemical characteristics 
of the propellant system. The method of calculation employed in the 

O 

present work was developed by Gordon and coworkers 0 , and adapted for use 
on the DRI Burroughs 5500 computer system. The reported calculations are 
based on 1967 JANAF Thermochemical Data.^ 

2.2 Calculation of Relative Mole Fractions . 

Propellant combinations with nitrogen tetroxide as an oxidizer and 
hydrazine, UDMH, or MMH as fuel can be classified as either H-N-0 or C-H-N-0 
atomic systems. Either classification implies that all combinations of 
compounds may be present that include these atoms and for which thermochem- 
ical data are available in the program. For the H-N-0 system the list of 
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possible reaction products is relatively simple, since the major products 
for this system include water, OH, NOg and NO. 

For the C-H-N-0 system the list of possible products is much more 
extensive. Typical calculational results for the propellant combinations 
involving hydrazine, UDMH, MMH and corresponding to the H-N-0 and C-H-N-0 
systjems respectively are presented in Figures lto l4. For the propellant 
combinations considered, about 30 individual species are possible. Of 
this number, only 23 are present within the defined minimum relative mole 
limit of 0.000005 moles. Calculations have been performed using a range 
of 0/F mixture ratios from 0.10 to 10.0 over a range of pressures from 
10,000 down to 0.3 psia. 

The pressure limits were chosen so that they would encompass the 
pressure regime between the initial start pressure and the high pressure 
transients caused by rapid shock heating during a gas phase detonation. 

The large variation in o/F was necessary in order to allow the estimation 
of the effect of variations of mixture ratio caused by flash vaporization 
or freezing lag effects. 

Although the major species could be estimated from the available 
literature, because of higher spectroscopic transition probabilities for 
some of the intermediate chemical species spectroscopic observations of 
the minor species was anticipated. The thermodynamic calculations outlined 
in Figures lto l4 indicate that a large number of intermediate species may 
be present under fuel-rich conditions. . Such conditions would be expected 
to exist in the vicinity of the fuel droplet. For an oxidizer-rich com- 
bustion system the calculated distribution of combustion products is much 
less complex. Significant amounts of NO 2 would be predicted in the vicinity 
of an oxidizer droplet. 


10 



PRESSURE = 0.3 PSIA 



1196 2472 2469 2353 2216 1908 1402 

TEMPERATURE, °K 


Figure 1 


11 


N 2 H4 / N TO 

PRESSURE = 6PSIA 


10 



0 

1.0 

2 0 

3.0 

4.0 

5.0 

0/F 

6.0 

70 

8.0 

9.0 

10 

1196 

2718 

2715 

2524 

2318 


1928 




1402 


TEMPERATURE, °K 

Figure 2 


12 


n 2 h 4 /nto 

PRESSURE = IOO PSIA 



1197 2929 2936 2645 2374 1936 1402 

TEMPERATURE, °K 


Figure 3 


13 


N 2 H 4 / NTO 
PRESSURE = IOOO PSIA 


10 



0 

1.0 

2.0 

3.0 

4.0 

5.0 

0/F 

6.0 

7.0 

8.0 

9.0 

10 

1203 

3053 

3083 

2708 

2399 


1940 




1402 


TEMPERATURE, °K 


Figure 4 


14 


n 2 h 4 / nto 

PRESSURE = 10,000 PSIA 


10 



. NH 

io~ 5 Li ■ i i i i i l 


0 

1.0 

2.0 

3.0 

4.0 

5 0 
0 /F 

6.0 

7.0 

8.0 

9.0 

10 

1244 

3122 

3183 

2745 

2412 


1941 




1402 


TEMPERATURE, °K 


Figure 5 


15 


PRESSURE = 0.3 PS I A 


10 


- H 



0 

1.0 

2 0 

30 

40 

5.0 

60 

70 

80 

9.0 

10 






0/F 






865 

2209 

2549 

2524 

2961 


2309 




1923 


TEMPERATURE, °K 


Fi gure 6 


16 


PRESSURE = 6.0PSIA 


i o 


- H. 



0 

1.0 

2.0 

30 

4.0 

50 

60 

70 

80 

90 

10 






0/F 






952 

2307 

2840 

2801 

2701 


2462 




1945 


TEMPERATURE, °K 


Figure 7 


17 


PRESSURE = 100 PSIA 



TEMPERATURE, °K 


Figure 8 


18 


PRESSURE = 1000 PSIA 



1.0 


2 0 


30 


40 


50 


60 


70 


0 0 


90 


10 


1259 2362 


0/F 

3363 3303 3079 2623 

TEMPERATURE, °K 


1957 


Figure 9 


19 


PRESSURE = 10,000 PSIA 


10 



Figure 10 


20 


UDMH/NTO 
PRESSURE = 0.3 PSIA 



814 1908 2537 2561 2519 2403 2111 

TEMPERATURE, °K 


Figure 11 


21 


UDMH/NTO 
PRESSURE = 14.7 PSIA 


10 



Figure 12 


22 


UDMH/NTO 
PRESSURE = IOO PSIA 


10 



Figure 13 


23 


UDMH / NTO 
PRESSURE = IOOO PSIA 



1225 1934 


3238 3426 3283 


2897 


2220 


TEMPERATURE, °K 


Figure 


14 


24 







SPECIFIC IMPULSE ~ FKOtBH 





III. MECHANISMS AND RATES OF CHEMICAL REACTIONS IN THE LIQUID 

AND GASEOUS PHASES . 

3.1 Introduction 

The thermodynamic analysis outlined in Section II provides a power- 
ful tool for predicting the products formed during the ignition and com- 
bustion processes under equilibrium conditions. In the ignition of a 
pulse mode engine, kinetic factors such as the rate of injection, flash 
vaporization, mixing and chemical reaction rates strongly modify such 
thermodynamic analyses. 

It is useful to review the reported chemical reaction mechanisms 
involving the decomposition of nitrogen- tetroxide and hydrazine, since many 
of the expected phenomena are attributable to the formation of metastable 
chemical intermediates which would not be formed if equilibrium conditions 
were present. These reaction mechanisms have been studied by a number of 
individual investigators One of the more comprehensive summaries 

of the rates of individual reactions associated with the oxidation of 
^2^4 ^ NOg "*- s ^ Sawyer'*--'-. In this summary, previous work in the 

field is reviewed, as well as data from a stirred flow reactor developed 
at Princeton x . The Princeton work can be interpreted that the reduction 
of NO 2 occurs as a two-step process in which the first step is a reduction 
of NO 2 to NO. The overall reaction mechanism suggested by Sawyer, based 
on experimental results from the stirred flow reactor, indicated a first 
order dependence with respect to the concentrations of both NgH^ and RO 2 
for the formation of NO. NO remains as a stable product for a relatively 
long period of time and subsequently reacts further to form N^ and H^O. 



Since the reaction delay for this second step is of the same order as the 
entire pulse-mode engine firing period, emphasis in the present review is 
placed on the first stage reduction step of this reaction sequence. 

3.2 Liquid Phase Reactions 

Liquid phase reactions between and have also been in- 
vestigated . Because of the extreme reactivity of this chemical 

system, the development of liquid phase study methods has been exceed- 
ingly difficult. Therefore, most of the recent investigators in the field 
of propulsion have assumed (rather than identified) the reaction products 
and have emphasized the effects of heat addition on rates of evaporation of 
liquid droplets. Because the currently assumed reaction mechanisms used to 
describe the thermal effects are somewhat doubtful, emphasis in this section 
is placed on a review of important liquid phase chemical reactions. Of 
significance is the Curtius reaction 

+'hno 2 — *hn 3 + H + + h 2 o 

which suggests that hydrogen azide would be produced. Also, it is noted 
that the oxidation of hydrazine by an oxidizing agent containing nitrogen 
usually produces HU .' Therefore, this review indicates that HN^, 
and HgH^-NO^ would be expected products of reaction for liquid phase re- 
actions occurring either in very dilute reaction mixtures, or at low 
reaction temperatures. 

3 

Perlee and coworkers , have identified hydazinium nitrate as the 
primary explosive present in engine residues using infrared spectroscopy. 



A verbal report from an in-house program at the Aerospace Corporation 
suggests that hydrazinium nitrite may also be present. From the point of 
view of gas phase reaction kinetics, hydrazinium nitrite would be expected 
to form in a direct bimolecular reaction between a collisionally stabilized 
HONO and a hydrazine molecule. This reaction would be expected to be 
reasonably efficient' and to produce the NgH^NOg intermediate. The forma- 
tion of the hydrazinium nitrate in the gas phase is more difficult to 
explain since it would require the production of HNO^ some point in the 

kinetic reaction chain. A possible mechanism for the formation of HNO^ in 
the gas phase would be the reaction 

HONO + N0 2 — > H0N0 2 + NO 

Although intermediate compounds of this type could form via a gas 
phase reaction mechanism and would subsequently condense, an alternative 
possibility is that hydrazinium nitrate forms in the liquid phase reaction 
between a major combustion product such as water, and condensed NgO^. Since 
an azide similar to hydrazoic acid was also found as a volatile reaction 
by Dauermann, et al, during a study of liquid phase reactions, it is neces- 
sary to develop a kinetic mechanism which can allow this formation to occur 
as a natural consequence of other reactions. A reaction scheme suggested 
as being reasonable involves the following sequence of reactions: 


H O + N 0, HNO-j + HNO (3-2-1) 

2 2 4 3 3 

H 2 H 4 + HNO 3 -♦> NgH^-NO^ (3-2-2) 

H H NO + HNO -+ HN + HN0-, + 2H 0 (3-2-3) 

^ P 3 ^ 3 < - 

+ N 2 H 4^ N 2 h / 3 (3-2-4) 
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In this reaction sequence, (3*2-3) is the Curt ius reaction. It is noted 

that two highly unstable compounds, HN, and ]OOL,, are predicted. Direct 

b 2 5 3 

chemical evidence for these two intermediate species. is described in detail 
in a later section. One interesting feature of this mechanism is that two 
moles of hydrazinium nitrate would be formed for each mole of water or 

consumed, which would suggest the formation of relatively large amounts 
of this intermediate product. 

An alternative direct reaction between hydrazine and nitrogen tetrox- 

ide , 


W + NHjNO + N g 0 + N 2 + H 2 0 (3.2-5) 

may also occur. Support for this latter reaction is from a wet chemical 
analysis performed, by Dauerman and infrared data obtained recently at the 
DFVRL in Germany. Since the NH^ infrared peak is masked by other reaction 
products there is some doubt as to the product identification. 

Hydrogen azide could disappear via Thiele's reaction: 


NH + HNO — ► N +NO+HO 
3 2 2 2 2 


(3.2-6) 


However, this requires a relatively high hydrogen ion concentration which 
is unlikely in this case because of the high percentage of unreacted hydrazine 
in the residue materials. 

3-3 Mechanisms and Rates of Gas Phase Chemical Reactions 

The kinetic mechanism proposed by Glassman and Sawyer 1 ^ f or the 

reaction between N^H^ and NC> 2 involves a two-stage reduction of NCk,. The 

first step of this reduction is a reduction of N0 2 to NO which occurs at a 

• . 1 

fast rate. This is then followed by the subsequent further reduction of NO 



to form N and water. The stage II reduction occurs at a slow rate and in 
2 

the flow reactor at Princeton was well separated from the stage I reduction. 
Since the ignition and combustion of a pulse-mode rocket engine will 
typically use a total hot firing cycle between 10 to 20 milliseconds the 
emphasis in this analysis will be placed on reaction involving stage I. It 
is assumed for this analysis that, both because of the short duration and 
because of the rapid quenching which occurs under space simulated rocket 
testing, the stage II reaction does not have time to proceed to completion. 
The Glassman-Sawyer mechanism depends primarily upon the stripping 


of hydrogen atoms from N^H^ and molecules. Studies by Gray and Thyne 

of H-atom abstraction from N-H bond by CH^ radicals demonstrate the sus^ 
ceptibility of hydrazine to attack. They found hydrazine to be 300 times 
more reactive than ammonia in this regard. Therefore, Glassman and Sawyer 
assumed that the initial step involved the abstraction of hydrogen from 
hydrazine to form H0N0 which subsequently decomposes to form OH and NO. 

A series of further abstraction was then followed, as outlined in the fol- 
lowing reaction sequence: 
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kl 

NgH^ +. N0 2 -*« N 2 H 3 + H0N0 
k2 

H0N0 — *> OH + NO 

k3 

N 2 H + N0 2 — *N 2 H 2 + H0N0 

k4 

N 2 H^ •+ OH — *N 2 H 3 + H 2 0 . 

N 2 H 3 + OH -^-N^ + H 2 0 

Additional hydrogen abstraction reaction by atomic hydrogen 


n 2 h 4 + .h 

N 2 H 3 + H 


n 2 h 3 + h 2 


N 2 h 2 


+ h„ 


(3*3-1) 

(3*3-2) 

(3*3-3) 

(3*3-4) 

(3*3-5) 

(3*3-6) 

(3*3-7) 


31 



and termination reactions 


. N 2 H 2~* N 2 + H 2 
H 2 + OH -* E 2 0 + H 2 

N0 2 + H-— ►OH + NO 

were other reactions also postulated by these authors. 


(3-3-8) 

(3.3-9) 

(3.3-10) 


Taking a simplified model based on the reactions (3*3-1) - (3-3-5) 
are assumed to be rate controlling, and then making the usual steady state 
assumption which is required by the principle of detailed balancing, the 
concentrations of the intermediates -are found to be given by the following 


expressions: 


0 ' 2 h 3 ] = (§;) 5 t " 2 H 4 ] 


[hono] = 


+ ^ k l k 3 k 4^2 


[H° 2 ] [n 2 h 4 ] 


(3.3-11) 


(3.3-12) 


(3.3-13) 


The reaction rate for hydrazine is then 


[n 2 h u ] 


= - k l + 


( 2)1 ^ ^ 


(3.3-1^) 


which has a first order dependence on both hydrazine and nitrogen dioxide. 
This first order dependence was observed in the stirred flow reactor at 
Princeton, which confirmed the proposed mechanism. 

A rate expression for N0 2 comsumption can be developed 

d [nOo] k l + ■ ( k l k 3 k 4) 2 ( 

— — = ‘ r [ H 2 H J (3-3-15) 
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which exhibts a first order dependence with respect to the initial re- 
actants. Again the type of dependence was observed experimentally in the 
stirred flow reactor and tended to confirm the mechanisms. Similarly the 


rate expression for nitric oxide buildup is given by 

. \ 

d [no] , \ ( kik^k^ 


k 2 [HONO] = 


dt 


k l + 


[NO ] [ N 2 H 4 ] (3.3-16) 


which is again first order with respect to both reactants. From available 
rate data this formulation appears to hold for high temperature, high pres- 
sure reaction systems. 

In the above scheme it is interesting to note that no account is 


taken of the termination reaction N H ■ 

2 2 


N + H which Glassman and Sawyer 
2 2 


suggest. If we assume that this dissociation reaction does not proceed to 
completion but that rather a condensation occurs 


xN 2 H 2 


°yy. 


(3.3-17) 


then (N H ) would appear in the condensed phase as a stable product. The 
2 x 

kinetic expression for its formation is 

• / • 

d [KjHj] _ 


dt 

d [ M a H a] 
dt 


= k 3 [n 2 h 3 ] [no 2 ] + k 5 [n 2 h 3 ] [oh] 

(k 3 k x k4) \ 


± k. 


[%HU] [BO,] 


(3.3-18) 


(3.3-19) 


•which is also first order with respect to N0 2 arid N^H^. This expression 
is identical to that for the consumption of NOg, indicating that the rate of 
formation of an explosive intermediate can also be correlated with the first 
order consumption of either of the reactants.. This buildup would be expected 
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to occur more readily at low temperatures, since the termination reaction 
( 3 - 3 - 8 ) must involve either a four-center decomposition or a hydrogen trans- 
fer followed by a three-center reaction. Data from the stirred flow re- 
actor cannot discriminate between a termination and condensation reaction 
by following the reactant decay rates unless a plot of temperature vs. 
energy released is also available. 

Low temperature conditions exist during the vacuum ignition of a 

rocket motor, and therefore it is predicted that substantial N^H^ could be 

formed during the preignition phase. Polymerization of N H 0 to form (NH) , 

2d x 

(where x > 2), produces triazene and tetrazene derivatives. Therefore, it 
is predicted that such compounds could be formed by gas phase preignition 
reactions. Evidence for such intermediates has been found in the liquid 
residue from rocket engines during this program. Data from the stirred 
flow reactor cannot discriminate as to whether this polymerization is occur- 
ring. 

At low pressure, reaction (3.3-2) is favored. However, in the re- 
actions where local pressures may be high, collisional stabilization to 
form nitrous acid would be expected. 

HONO* — *- HONO (3.3-20) 

If we assume that nitrous acid is collisionally stabilized, a simplified 
Glassman-Sawyer reaction mechanism results. 

(3-3-1) 
(3.3-3) 
(3.3-4) 
(3-3-5) 


N 2 H 4 

+ NO 

2 

—t N 2 H 3 

+ HONO 

n 2 h 3 

+ no 2 

-* N 2 H 2 

+ HONO 

N.H, 

+ OH 

— *-N H 

+ HO 

2 4 


...23 

2 

N H 

+ OH 


+ H 0 

2 3 


2 2 

2 
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In this simplified mechanism HONO'and are considered as stable end 

products formed during the preignition period of operation. Using the 
steady state assumption 


d [n 2 H 3 ] 

= 0 = k x [ N 2 H 4 ] [N0 2 ]- k 2 [N 2 H 3 ] [N0 2 ] 

d."fc 


(3.3-21) 


the rate of reactant consumption becomes 


d ["A] 

dt 

d [N0 2 ] 
dt 


= - kn 


[N 2 H 4 ] [no 2 ] 


2k L [W 2 H 4 ] [N0 2 ] 


(3.3-22) 


(3.3-23) 


for the fuel and the oxidizer respectively. This simplified mechanism 
can be extended to provide a basis for the formation of NO and E^O by two 
consecutive hydrogen abstractions by NO^ involving either ( 3 . 3 “l) or 
( 3 « 3 “ 3 ) followed by 


HONO + N H, ^ H 0N0 + N H (3. 3-24) 

2 4 2 2 3 

hono + n 2 h^ !si*. h 2 ono + n 2 h 2 (3.3-25) 

h 2 ono h 2 o + no. (3.3-26) 

An analysis of this scheme shows the same first order dependence that was 
present in the G-S mechanism. 

A third mechanism involves the stabilization of HNO intermediates 

by the following series of reactions 

W 2 H 4 + N0 2 — * HONO + N g H (3.3-1) 

HONO — ► NO + OH (3.3-2) 

NO + N H, i— ♦> HNO* + N H (3-3-27) 

2 4 2 3 
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with the HNO* becoming stabilized. Evidence for this type of stabilization 

is the observation of (HNO) emission spectra and the presence of (HNO)^ ' 

infrared absorption bands in the condensed residues. Although this radical 

is unstable it has an upper state electronic "lifetime" of about 0.1 second 

which is an order of magnitude longer than the pulse-mode combustion process 

For this reason it would be possible to form substantial quantities of 

(HNO) under the reaction conditions present in the pulsed rocket engine, 
x 

As described later, during the high altitude simulation tests the 
nitrogen tetroxide forms a fine mist whereas the fuel stream remains as an 
integral liquid jet. This process results in a high concentration of 
NOgCg) surrounding the fuel droplets. As the fuel droplets evaporate, a 
stepwise hydrogen abstraction can occur resulting in the possibility of 
completely stripping hydrogen from the hydrazine. 

n 2 h 2 + no 2 H0N0 + iyi (3-3-28) 

N H + NO Np + H0N0 (3*3-29) 

N 2 H + N 2 H 2 ^ N U H (3.3-30) 

This suggests that an oxidizer lead could lead to the production of complex 
gas phase azide compounds. Experimental evidence for such compounds was 
found in chemical analyses of residues in this program. 

3.4 Survey of Specific Rate Constants Available from Literature 

A literature search was performed in order to document the range of 
specific rate constants that have been investigated for the gas phase re- 
actions involving nitrogen tetroxide and hydrazine. Because of increased 
complexity and limitations in time a similar search was not performed for 
MDH and UDMH decomposition reaction. The results of this effort are con- 
tained in Appendix A for nitrogen . tetroxide; and hydrazine and their de- 
composition products. 
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NO z REDUCTION (1st Stage) 

N 2 H 4 + NO z — N 2 H 3 + HONO 

HONO OH + NO 

N 2 H 3 + NO z N 2 H 2 + HONO 

N 2 H 4 + OH — *■ N 2 H 3 + H z O 

N 2 H 3 + OH *»" N 2 H 2 + H z O 

Postulate 

NO + N 2 H 4 — HNO + N 2 H 3 

HNO (HNO) x 

NO + N 2 H 4 NOH + N 2 H 3 

NOH H 2 N 2 0 2 


Figure 17 ^2 Reduction .Sequence 


HYDRAZINE DECOMPOSITION 
M 

N 2 H 4 NH 2 + NH 2 

N 2 H 4 + NH 2 — — N 2 H 3 + NH 3 

N 2 H 3 NH 2 + NH 

N 2 H 3 + NH 2 N 2 H 2 + NH 3 

n 2 h 2 — - n 2 + h 2 

Postulate 

N 2 H 2 (NH) X (x > 2) 


... - * m y * - 7 

Figure;' 18. '■■Hydrazine, Decomposition Sequence 


hydrazine/nitrogen tetroxide 

(Stage I) 
k 

N 2 H 4 + .N0 2 — U- N 2 H 3 + HONO 

HONO HO + NO 

N0 2 + N 2 H 3 HONO + N 2 H 2 

N0 2 + N 2 H 2 -^ i — N 2 H + HONO 
N0 2 + N 2 H — Nj + HONO 

N 2 H 3 NH + NH 2 

NH + N 2 H 4 NH 2 + N 2 H 3 

NH + NO z HONO •+ N 

N + N0 2 - 1 *- NO + NO 
k 

NO + N0 2 — N 2 0 3 NO + NO z 

NO + N 2 H 4 — HNO + N 2 H 3 

k, 2 

NO + N 2 H 3 — ^ HNO + N 2 H 2 
NO + NH 2 - k - 13 * - HNO + NH 
NO + NH * • HNO + N 

NH 2 + N 2 H 4 — ^ N 2 H 3 + NH 3 
NH 2 + NO z HONO + NH 

OH + N 2 H 4 - kl7 > N 2 H 3 + HOH 
Overall 

N 2 H 4 + N0 2 - HNO + H z O + N 2 

Figure' 1'9. Overall Hydrazine/Nitrogen 
• \ . . ^-Tetroxide ' React ion Sequence . 

)>o 




REACTION RATE 
(degrees Kelvin) 
a 300 
b 1000 
c 2000 


Figure 20 


Morphological Kinetic Rate Diagram for the Decomposition of Hydrazine 





IV. 


EXPERIMENTAL PROCEDURES 


4.1 Introduction 

An experimental test program was developed in order to investigate 
the physical and chemical processes responsible for the observed detonations 
in pulse mode engines. The purpose of the test program was to study the 
buildup of residue material in rockets under laboratory conditions. The 
major design criterion for this work was an operational thrust chamber 
pressure of 100 psia and short pulse capability in order to simulate high 
altitude operation of the Marquardt RL-10 reaction control system rocket 
engine. 

To perform these tests, a highly instrumented 2 dimensional rocket 
engine and a high altitude simulation test facility were designed and built. 
The experimental test chamber had a- volume of about 550 cubic feet and an 
altitude simulation capability of 200,000 ft. (0.15 Torr). 

In order to allow a direct comparison of data obtained in this pro- 
gram with data reported by a concurrent test program at the Bureau of Mines, 
the thrust chambe-r and injector configurations were similar to those used 
at the Bureau of Mines test facility. Therefore, the pulse mode rocket 
developed for use on this program was a two-dimensional thrust chamber 
utilizing a one-on-one injector. Development of the thrust chamber also 
involved the design and construction of relatively complex propellant 
handling, and propellant flow control systems. Instrumentation for this 
test facility required a, very high response tape recording system in order 
to record high frequency transducer data and to monitor fast response 
spectrophotometric output. 



4.2 Test Facilitie 


4.2.1 Two Dimensional Rocket Engine 

A configuration similar to the Bureau of Mines' engine was used for 
these studies. Thrust chamber geometry, injector impingement angles, 
impingement distances, and inlet orifice dimensions were identical, in 
order to allow comparison of data generated between the two concurrent 
programs . r However, a number of differences exist in the configuration of 
this engine and the Marquardt R4-D. The major difference involves place- 
ment of a single one-on-one doublet injector along the center axis rather 
than a set of eight doublets arranged on a concentric circle normal to 
this axis. This arrangement reduces the effective mixing and allows the 
propellants to impinge on the aft chamber walls. A second differen.ee is 
that fuel weepers used for thrust chamber cooling in the Marquardt engine 
were not incorporated in the design and impingement of the fuel onto cold 
chamber surfaces near the injector plate did not occur. The Reynolds 
number of the propellant flows is calculated to be about 2.8 times greater 
than that of the operational engine although, since both should be in a 
stable turbulent flow regime, no significant effect is expected from this 
factor. 

The main elements of the design departure from the Bureau of Mines 
engine can be summarized as follows: 

1. A 6061 aluminum thrust chamber fitted with viewing ports 
replaced an all plexiglass engine. 

2. Direct electronic monitoring of reactant flow rates replaced 
photographic flow monitoring system. 



3 . A valve system which enable recycling under vacuum conditions 
in addition to precisely controlled reactant injection re- 
placed a burst diaphragm flow controller • 

4. Thrust chamber structures capable of withstanding internal 
detonations were employed. 

The aluminum housing was chosen because the number of firings required a 
rugged, nonreactive, easily machinable, construction material. The window 
orientation was chosen to maximize the optical depth of the reaction chamber 
in order to increase the emission and absorption signals. Flow rates were 
monitored using a linearized electronic readout on a high-speed tape system. 
To allow this readout, the piston in each reactant reservoir was fitted with 
a core of permeable material that affected the inductance of a nonuniform 
coil wound along the wall of the reservoir. The inductance of the coil was 
then a function of the core position and could be related to the instantan- 
eous flow rate of the reactants. 

Design drawings of some of the major components are presented in 
Figures 21 through 25 . Figure 21 shows the thrust chamber. This is a 
two-piece sandwich construction in which the individual components are 
bolted together. Window ports are located to view the impingement point of 
the fuel and oxidizer streams. The materials of construction which were 
used are a 6o6l aluminum alloy for the thrust chamber and quartz optical 
windows. Figure 22 shows the injector block. Impingement angles of 30° 
and 40° are used for fuel and oxidizer propellant streams. . Figure 23 shows 
the valve assembly and fuel inlet system. Valve action is performed by a 
rapid acceleration of the Teflon piston via a rapid pressure release on one 
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Figure 2-U» . Schematic Diagram of Propellant Handling System and Experimental 
Test Chamber 




Detail of the Construction of the Reservoir Monitor (Nonlinear 
Inductor) and Piston Design 



Time, I mil lisecond /division 

Example of Transient Response Obtained From Reservoir 

Volume Monitor 

Figure 25. Reservoir Monitor and Piston 
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end of the piston. The valve is essentially a piston that contains a 
mating passage that matches the injector passage when the piston is trans- 
lated into position. The Teflon piston is actuated pneumatically in order 
to minimize piston mass and to obtain an adjustment actuating force. The 
system features ppositive sealing which allows prepressurization of the 
reactant reservoirs, and the filling method prevents the formation of a 
vapor pocket within the reservoir. The design also isolates the reactant 
supplies from the chamber during the filling cycle and from both the 
chamber and reservoir during firing.' This approach allows the motor to 
be recycled while under vacuum conditions together with the advantage of 
the fixed-volume injection scheme. The reactant reservoirs are shown in 
Figure 25. 

4.2.2. Vacuum Test Facility . 

A vacuum test facility was assembled and evacuated to a pressure 
corresponding to an altitude of 200,000 feet. A schematic diagram of 
this facility is given in Figure 26. A stokes 212G (130 cfm) and a 
Kinney KMB-230-130 pumping system was used. A detailed view of the optical 
port for this facility is given in Figure 27. This facility consisted of a 
large vacuum chamber (490 cubic feet) and a thick-walled test section 
(20 cubic feet) connected by a l6-in. flange. The test chamber has 
\ in. steel walls which are capable of sustaining moderate blast waves. 

This chamber was fitted with three 8- in. observation ports . The large 
chamber is of riveted steel construction and is used to provide a ballast 
volume for the test chamber ... 



- Ballast Chamber 


Test Rocket Entry Port (12 FI.) 

Observation Ports (8” FI.) 

Power/water Service Station 



-Kinney Vacuum Pumps (1st & 2nd Stage) 


Stokes (MicroVac) Supplementary Pump 


'Vacuum Gage 


Figure 26. Sketch of Rocket Test Chambers and Vacuum System 
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4.2. 3 Propellant Handling System 

Schematic diagrams of the propellant handling system are shown in 
Figures 29 through 31 • Figure 31 also shows the relative placement- of 
engine hardware within the test chamber. To achieve fast dynamic response 
for. valve actuation, helium gas is used as a pressurant. The main fuel 
and oxidizer handling systems are mounted directly on a demountable vacuum 
test chamber face plate. ~ 

4.3 Instrumentation 

- A variety of techniques were used to monitor the chemical and 
physical processes that occurred during the pulse-mode operation of the 2-D 
rocket engine. Pressure measurements were made using two Kisler 601A 
pressure transducers. These transducers were located in the upper wall of 
the thrust chamber. One transducer was emplaced near the point of propel- 
lant impingement; the second transducer was located in the aft end of the 
chamber just before the converging throat section. The output of these 
transducers was simultaneously displayed on a Tektronic 547 or 555 oscil- 
loscope and recorded on a l4 channel Mincom. M34 high-speed tape recorder. 

To amplify the transducer signal two charge amplifiers were used. The 
frequency response of the tape recorder was dependent on whether direct or 
FM mode of operation was used with a frequency response of 80 kilohertz and 
600 kilohertz being the maximum available from the FM and Direct mode 
respectively. In order to reduce a zero point drift of the transducers 
caused by thermal shock, a very thin film of grease was applied to the face 
of each transistor prior to hot testing. 



v_n 

VJI 



Fuel € Oz.«<li ?er 
Resevotr Le>o.<t»r\Q 

■j Rsfcon Vo.ft/e. 

fVess o- r i 0.4. i on 

. 

~v Piston Ci-ee.4) I'kiu’e, 

I /^ttu/aL-trion 



Copper < Br<uss 




Figure 3 0. Schematic Diagram of Fuel and Oxidizer Pressurization Systems 
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Figure 31. Injector Block with Chamber Detached 


Propellant flow measurements are made using the specified nonlinear 
inductor described in the previous section. The frequency response of the 
inductor was limited by the frequency response of a set of preamplifiers to 
6 kilohertz. 

Spectral data were recorded using a number of different • spectrographs 
and spectrophotometers. Time integrated spectra in the visible region 
between 260C$ and 750C>2 are taken using a Bausch and Lomb 1.5 meter grating 
spectrograph. This spectrograph has a resolution of l4.763^/mm. Time 
resolved spectrographic data were obtained using a CINE quartz prism spectro- 
graph. This spectrograph is a specialized spectrograph that had been developed 
to record time resolved spectral data emanating from large scale detonations, 
and has a spectral response between 2200 X and 65 OOA, covering the ultraviolet 
and visible regions. However, because of its inherent low resolution only 
a limited amount of work could be performed with this spectrograph. A one 
meter Jarrell-Ash spectrograph was modified for use on a polychrometer . • 

This instrument uses a grating -of 1158 lines/mm with useful band width of 
o 

about oOOA. This instrument was modified to obtain time resolved short-band 
pass spectral records by incorporating a rotating quartz refraction prism 
in the optical path. A Warner Swasy Model 501 high-speed scanning spectro- 
photometer was used in the visible and infrared spectral region. 

A combination of emission and absorption spectra was used to follow 
the reaction processes occurring during the operation of the rocket thrust 
chamber. A controllable pulsed backlight source lamp was developed to aid 
in this analysis. A Hewlett Packard 213 pulse generator with a Xenon Arc 
Lamp which was collimated into a probing beam approximately 3/8 inch in 
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diameter was also designed and fabricated. In order to avoid overdriving 
the lamp for other than brief time periods a short duration, repetitive 
pulse was imposed on a continuously operated lamp. 

A typical oscilloscope output of the light system is shown in 
Figure 32a, together with a plot of the relative radiant output character- 
istics of the lamp as a function of wavelength. Typically, light pulses 
with a 15 microsecond half -width at a repetition rate of 5 kH were obtained. 
This repetition rate could be substantially increased. The light was run 
continuously at a low level and is pulsed to a transient intensity of 
approximately four times the steady state output. The intensity shift is 
shown in the oscilloscope output tracing. A line spectrographic record of 
this lamp in operation is shown in Figure 32b. 

High-speed photographic studies were performed using either a HyCam 
or an Eastman l6mm camera. Both cameras were used to obtain data at about 
3000 frames-per-second although the HyCam had a much higher framing rate 
capability. 

4.4 Cold Flow Testing 

Cold flow tests were performed to test the design of the overall 
system. These tests had three objectives: (l) to test the integrated 

control and propellant handling systems; (2) to aid in the development of 
operational test procedures; and ( 3 ) to test the operating characteristics 
and reliability of the propellant valving and control systems. As outlined 
previously, specialized valves were designed and constructed in order to 
allow recycling of the engine system without breaking vacuum, and to retain 
the concept of using fixed-volume piston-driven propellant reservoirs. 
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Figure 32. Performance of Pulsed Xenon Lamp Measured at 500 nrn 
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Precise and accurate propellant control in terms of lead, lag, and 
rate of flow were achieved with this valve system using water as a propel- 
lant simulant*. No leakage of propellant around the moving seals was ob- 
served under either ambient or high altitude pressure conditions. Figure 33 
shows a stroboscopic view of a typical cold flow test. Total propellant 
flow for this test is 2.5 cc. and total firing duration is 20 milliseconds 
with a 5 millisecond fuel lead. 

A relatively large number of cold flow tests (approximately fifty) 
were performed in order to determine the reproducibility of propellant 
handling systems and to outline necessary operational pressures required for 
specific flow rates.. Bqsed on these tests, it was determined that the fuel 
(or oxidizer) leads could be sequenced with high relative precision. Since 
the unique flow monitoring methods developed could be directly correlated 
with actual propellant flow a relatively precise knowledge of the exact 
propellant flow was provided. These tests indicated that the test design 
concepts, as developed," in fact provided a highly flexible and precisely 
controlled propellant handling system. 

4.5 Single Propellant Tests 

Tests using single streams of individual propellants were performed 
in order to obtain background data concerned with the spray patterns and 
spray droplet size distribution functions. Results of these tests indicated 
that the initial spray patterns were strongly dependent on propellant 
physical properties. 

During cold flow testing two problems developed with the recommended 
grease used to seal the moving "o" ring on the propellant reservoir. The 


61 




Figure 33. Strobe View of Water on Water Injector 
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first problem involved the lack of adherence of the grease to the piston 
wall surfaces, caused by an apparent setting of the grease in a semi- 
regular microcrystalline structure. The second problem involved the 
chemical attack of the liquid oxidizer on the grease. A series of com- 
patibility tests was therefore initiated and a change of lubricant was 
implemented. The Dri-Lube* grease which undergoes only limited chemical 
degradation upon contact with liquid WTO appears to be a satisfactory 
lubricant . 


*Dri-Lube ; Type 822; MMSN 306A, Lox, NTC, Hydrazine Compatible. 
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V. 


EXPERIMENTAL RESULTS 


5.1 Introduction 

• Using- the test facilities and instrumentation described above, 
three series of rocket tests were performed. The first series involved 
cold flow tests in which simulated propellants were used in order to 
verify the design and operation of the rocket engine and propellant handl- 
ing systems. During these tests calibration of the fuel flow rates and 
verification of electronic sybsystems were performed. 

A second series of tests involved studies of the effects of high 
altitude simulation on propellant injection hydraulics. For this pur- 
pose a large number, of tests were performed using single propellants. In 
addition, a few selective hot firings were performed in order to demonstrate 
the combined effects for typical propellant systems. These tests were used 
to provide information on spray distribution patterns and to investigate 
the effect of changes in propellant physical properties on droplet vaporiza- 
tion rates. 

The third series was a hot firing test program. During this series 
of tests, flow, pressure and spectral data were recorded. In some cases a 
chemical analysis of the residues formed during pulse mode operation was 
also made in order to verify and extend the results of other investigators. 

A more detailed description of these tests together with a documenta- 
tion of typical results is presented in' this section. 

5.2 Cold Flow Tests 

An extensive series of approximately 150 "cold flow" tests were 
carried out in order to verify the design and operation of the 2-D rocket 
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engine and the propellant handling systems. A relatively large effort was 
directed toward the dynamic calibration of the propellant flow monitoring 
system. 

Calibration curves were prepared in order to correlate the propel- 
lant flow output signals with actual propellant injection rates. The 
calibrations indicated a high degree of repeatability (-2 'jo) between con- 
secutive tests. Figure 3^-b shows a signal obtained from a typical test 
run. In this oscilloscope record, the piston travel displacement is dis- 
played in the ordinate with time on the abscissa. As indicated, the time 
base is 1 millisecond per major division. The slightly curved shape of the 
flow response curve in the central region is a result of variations in 
induction coil spacing which were hand wound, and does not indicate a 
fluctuation in propellant flows. In reducing the data it is necessary to 
subtract this variation in order to obtain a true flow, response curve. 

Figure 34a shows a typical spray pattern obtained from these tests. 

- On this test the thrust chamber had been removed in order to allow better 
observation of the development of the spray pattern. As can be seen, there 
is a relatively even distribution of material throughout in the spray area. 
The slight asymmetry of this pattern is caused by a difference in momentum 
vector of the two propellant streams which results from the difference in 
density between water and nitrogen tetroxide. A crude calculation of this 
momentum imbalance effect was made and it was concluded that the resultant 
spray distribution for an actual "hot" test would be uniformly distributed. 
Considering the one-on-one type of injector used, this experimental result 
was considered to be within the range of variation of design parameters. 
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Figure 34. Spray Pattern and Oscilloscope Trace of the Dynamic 

Response for Propellant Flow and Pressure Transducers 
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To provide detailed information with regard to the injection 
hydraulics and the propellant sequencing parameters, the development of 
sprays during cold flow tests was monitored using high speed motion 
pictures. These photographs showed in considerable detail the physical 
phenomena which occurred during propellant sequencing. A relatively 
complete series of tests was conducted using water to simulate both fuel 
and oxidizer. A portion of a typical sequence showing some of the details 
involved in the injection hydraulics is given in Figure 35 • 

At the onset of propellant flow, a small filament corresponding 
to leading edge of the lead propellant is observed. In most cases this 
edge has either an upward- or downward momentum vector corresponding to the 
opening of the oxidizer or fuel valve respectively. Under high vacuum con- 
ditions the leading edge is developed into a "blossom" in which the propel- 
lant stream is effectively broken up. The reason for this transient blossom 
effect is not presently understood, but it may involve either a pre-nuclea- 
tion of the propellant by the dissolution of pressurant gases or a breakup 
caused by the increased turbulence when injectors having large l/d ratios 
are used since during valve opening the orifice has a rapidly changing 
cross section. In the latter case, it is probable that the L/D as well as 
the injector characteristics are substantially changed by the valve actua- 
tion. The propellant stream assumes the form of a highly turbulent column 
immediately after the onset of main propellant flow. 

Upon termination of propellant flow these photographs demonstrate 
an extremely sharp fuel or oxidizer flow cut-off. This is in contrast to 
most propellant engine systems in which manifold effects are present with 
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Figure 35. 


Vacuum Injection Behavior for N 2 H 4 at 100, 000 Feet ~ 3, 000 


68 







The concommittant dribble volume . The propellant injector developed on this 
program is essentially a zero dribble volume system which eliminates the 
formation of hydrazinium nitrate residues formed from the residual propel- 
lant emptying from injector manifolds and propellant lines. 

5-3 High Altitude Propellant Injection Studies 

A series of tests was performed in order to outline the characteristic 
behavior of hydrazine/NTO propellants under simulated high altitude startup 
tests. Forty-three tests were run using both single propellant flows and 
hot firings. For a given injection momentum there was a strong altitude 
effect on propellant stream disintegration. In order to permit better 
observation of the behavior during the preignition phase under vacuum start- 
up conditions, the thrust chamber was removed. Characteristically, under 
high altitude simulation the fuel stream remained integral throughout the 
flow region observed. This was in marked contrast to the behavior of the 
oxidizer, which underwent rapid disintegration to form a very fine aerosol 
at the face of the injector. 

Three typical spray patterns are shown in Figure 36 . As can be 
observed, the oxidizer fans out in a wide spray, whereas, the fuel (WgH^) 
retains its jet integrity even at very low pressures. This behavior can 
be predicted on the basis of vaporization theory providing essentially 
empirical correlation parameters are used. For a hot firing the effect of 
the large differences in vaporization is graphically shown by the last 
photograph. The hydrazine consists of a liquid column immersed in a fine 
aerosol of NTO. This characteristic behavior does not occur during ambient 
pressure testing. All photographs were taken on a HyCam camera operated at 
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N 2 H 4 3 mm 


N 2 0 4 10mm N 2 H 4 /N 2 0 4 100mm 

Figure 36. High Altitude Spray Pattern 


300 PPS using 2484 Eastman Kodak Film (ASA 2000). The event is frontlighted 
using a bank of three, high intensity lights. 

5.4 Hot Firing Tests 

A large number of hot firing tests using hydrazine/NTO propellant was 
performed. The major objective of these tests was to obtain spectral data 
which would provide information regarding the chemical processes responsible 
for the catastrophic overpressures. Approximately 90 tests were run at low 
ambient pressures. A number of detonations were observed during the course 
of testing. In most cases those detonations were nondestructive in that the 
engine hardware remained intact, but the force of the detonation was oc- 
casionally sufficient to completely shatter the ^ inch quarts viewing ports. 

5.4.1 Pressure and Flow Response Data 

Using the instrumentation described previously, data on pressure, 
flow, and spectra were recorded. An oscilloscope trace of the dynamic 
response for propellant flow and pressure transducers is shown in Figure 37 . 
The time base for this test is 20 milliseconds per major division with the 
time axis running from left to right. Because this record is typical of 
most of those taken, a detailed explanation of the trace is given below. 

Pressure data were recorded on the upper two traces of the oscillo- 
scope and measured at two locations. No significant time variation between 
the two locations were observed during any of the testing performed. In 
Figure 37 a deflection on the vertical scale (upwards) indicates an increase 
in pressure and a standardized vertical scale deflection of 500 psia/cm was 
used throughout these studies. In this record a pressure spike of about 
1100 psia was recorded although pressure levels well above 2000 psia were 
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Figure 37. Dynamic Responses for Propellant Flow and 
Pressure Transducer 
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observed during these tests. There is a zero point shift which occurs 
i mme diately after the start of the pressure signal. This shift limits the 
accuracy of measurement to about 15 percent with a greater uncertainty in 
measurement at the higher pressure levels. The cause of this shift is not 
known, although it may be attributed either to a thermal shock or to a 
contact resistance effect. A thin layer of grease applied to the surface 
of the transducer in order to reduce thermal shock was only partially 
effective. A series of tests in which a photoflash bulb was placed approxi- 
mately 2 inch from the transducer surface verified the hypothesis that the 
baseline shift is in some way related to a thermal shock caused by the im- 
pingement of radiant energy onto the transducer surface. Conditions were 
adjusted so that there could be no pressure fluctuation during this test. 

The larger shift in the second trace possibly may be attributed to increased 
thermal shock since it is located in the aft end of the thrust chamber where 
local temperatures are somewhat higher. This zero point shift causes a 
large uncertainty in the absolute value of the recorded pressure level. As 
can be seen by inspection, the baseline returns asymtotically to the original 
value after about 200 milliseconds. 

Propellant flow rates are recorded on the lower two traces. Fuel flows 
are recorded on the upper trace and oxidizer flows on the lower trace. For 
precise analysis it was necessary to take into account a slight fluctuation 
in signal in the central portion of the propellant reservoir caused by the 
inductor coil windings. A careful analysis of these records shows a fuel 
lead of 3.5 milliseconds. Seven milliseconds later a sharp reversal in the 
flow of both the fuel and oxidizer occurs. This fluctuation in flow cor- 
responds exactly with the pressure transducer records indicating a high 
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pressure ignition spike. The fine structure of this record also indicates 
a series of three flow perturbations, each of -which is directly correlated 
with a pressure fluctuation in the thrust chamber. These results demon- • 
strate that the instrumentation, as developed, can monitor direct inter- 
actions between the ignition process and the injection hydraulics. These 
fluctuations occur at a frequency of about 200 cycles per second and could 
therefore be described as a low frequency chugging. 

An insight into the transient flow characteristics at the onset and 
shutoff of propellant flow can be obtained by close examination of the 
shape of the flow response curve. At the onset of propellant flow the flow 
• curve has a round shoulder which indicates a constantly increasing flow rate 
during valve actuation. At the termination of flow a very sharp discon- 
tinuity is observed, indicating an almost instantaneous cutoff of flow. 

This correlates with visual studies obtained using high, speed photographic 
techniques. 

5.4.2 Infrared Analysis of Engine Residues 

The severe ignition overpressures sporadically observed during the 
vacuum ignition of NTO/A-50 propellants are believed to result from a 
detonation of the residue material formed in the chamber. Therefore, a 
careful chemical analysis was performed on residue samples obtained from 
hot engine firings. Figure 38 shows the infrared absorption pattern obtained 
from a. typical MMH/NTO residue. 

In addition to hydrazinium nitrates several other species were detected. 
These included triazene > hydrazoic acid (HN^) , hyponitrous acid 

(NgH^), and monomethyl triazene (CH^^ ) . These species were present in 


74 


MICRONS 

3 4 5 6 



Figure 38 . Typical Infrared Absorption Pattern for NTO/MMH Residue 






small amounts and. their significance is believed to be associated with 
their possible role in increasing the explosive sensitivity of the base 
material. Triazene, monomethyltriazene and hydrazoic acid were identified 
in both mass spectrometric and infrared analysis. Hyponitrous acid was 
identified in the infrared analysis only. For these compounds all reported 
bands between 2.5 and 15 microns were either specifically present or, 
alternatively, were positioned at a frequency corresponding to a major 
absorption band of hydrazinium nitrate and/or one of its derivatives. The 
single exception to these patterns was a weak combination band + v ^ 
for NgOo which was not observed. 

Mass spectrometric data (outlined in the following section) indicated 
sizable M/e ratios at 43, 45, and 59 which are correlated with HN^, 
and respectively. Secondary breakdown peaks could not be used 

because of a lack of compound reference breakdown data for the mass spec- 
trometer used. These data were obtained from the volatile material evolved 
from a sample of the residue which had been frozen using liquid nitrogen 
and then evacuated to pull off dissolved gases. Therefore, the mass 
spectral data provides information about the character of the gaseous 
species which are formed by reactions taking place in the liquid phase of 
the residues. In general the mass spectral patterns were well correlated 
with infrared absorption data. 

The relevance of the presence of the observed metastable compounds to 
the problem of ignition overpressures and detonation is that they can be 
expected to act as sensitizing agents. It is probable that these compounds 
are formed in different stages of the reaction. For example, it is hypo- 
thesized that HNO is formed as a. stable intermediate during the gas phase 
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reaction of NO^ with hydrazine. Subsequently, this dimerizes to form 
hyponitrous acid. Since a well known synthetic route for the formation 
of HN^ is the reaction of HgHgOg and IT in the cold, it is probable that 
HN^ is formed in liquid phase reactions between the unburned propellant 
and quenched intermediate combustion products during the shutdown phase of 
engine operation. 

5.^.3 Mass Spectral Analysis of Engine Residues 

Additional information regarding the chemical composition of the 
engine residues was obtained by a mass spectral analysis of the volatile 
components emanating from the residue formed after a MMN/nTO engine firing. 
Sample material was obtained from the thrust chamber wall after a firing in 
which a severe detonation had occurred at the observation window port. 
Because of the mode of fracture at the window it is believed that the 
material collected represented a sensitized but undetonated liquid explo- 
sive. 

As shown by Per lee, an overpressure of the order of 22 kilobars is 
required to detonate hydrazinium nitrate in an "unconfined" space. No 
evidence was found from transducer data that pressure levels of this order 
were present in the chamber. The sample was held at low pressures for 
several hours at room temperature conditions. Although the sample was held 
under an inert atmosphere there is a possibility that air oxidation in 
addition to slow liquid phase reactions occurred. 

A typical mass spectral record is given in Figure 39 • This record 
plots the observed peak height vs. the mass-to-charge ratio. Corrections 
have been made for the background. Two scale sensitivities are given at the 
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abscissa. Of specific interest to this research are the peaks observed 
at M/e values of 43, 45, and 59 which are correlated with hydrazoic acid, 
triazene and monomethyltriazene respectively. The ratio of the 43/42 peak 
for pure HN^ reported by Ray et al, is 10:1 which is within the experimental 
range observed during these studies. No information was found for the 
spectral breakdown pattern of triazene and monomethyltriazene. A detailed 
analysis of the breakdown pattern at lower M/e values was not performed 
since the mixture has a complex breakdown pattern in the lower M/e region. 
Because the mass spectrometer used a magnetic scanning principle there are 
substantial differences in M/e peak sensitivities, compared to values given 
in the American Petroleum Institute Tables. 

From these analyses it is concluded that there is evidence for the 
presence of hydrogen azide, triazene, and monomethyltriazene in the liquid 
residues. These compounds are known to be relatively shock sensitive and 
may act as sensitizing agents, i.e., as explosive initiators. Because of 
the delay in sampling and analysis it is not known whether these are formed 
on a time scale similar to the engine pulse width or at relatively slow 
rates. It would be desirable to monitor the gases evolved from the residues 
using a time -of -flight mass spectrometer directly attached to the thrust 
chamber . 

5.4.4 Emission Spectra Obtained from Hot Firings 

Figure 40 shows a typical time integrated spectrum from an MMH/NTO 
pulse mode engine firing. The spectrum was obtained using a grating Bausch 
and Lomb spectrograph having a resolution of l4.753S/mm. Because of the 
brief duration of combustion, it was necessary to develop a high performance 
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Figure 40. 



film-developer combination. The best compromise between grain size and 
speed was found to be Eastman Kodak 2485 film using an m 62 developer. 

This combination provided an estimated ASA rating of about 1200. Figure 4l 
shows a densitometer tracing of the same spectrum. As can be determined 
from these figures, there is a relatively complex spectral emission pat- 
tern exhibited by the reaction. Atomic Na and Li are prominent and serve 
as a calibration check. The gross structure consists of a series of bands 
degrading to the red with band heads at approximately 4060, 4280, and 4650A, 
together with a series of weaker emission bands extending from about 3200 
to the upper end of the red cutoff limit of the film. 

The analysis of the spectrum is divided into two parts. The first 

part is concerned with the identification of the principal emission, bands , 

while the second part is concerned with the identification of the weaker 

fine structures. This analysis indicates that the species responsible 

for the three main band structures is an electronically excited metastable 

state of NO . The fine structure or the weaker bands, which can be detected 
2 

by using a densitometer, consists of a large number of individual species 
and for the most part represents transient metastable compounds. 

The principal emission structure in the region between 4000 and 
500 oX can be correlated with an N0^ fluorescence between a metastable upper 
electronic state and the ground state. One major problem in making this 
assignment is the unusual clarity of the band structure. This unusual 
clarity is difficult to explain unless it is assumed that the usual ab- 
sorption bands are limited to a specific transition associated with a 
chemiluminescence reaction. Previous studies in this wavelength region have 
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assigned the broad continuum in the region derived from the same fluores- 
cence, and commonly referred to as the afterglow, to a reaction product 

* 

formed, by the reaction. 0+N0 — ►NC^. More recent studies indicate that 

this continuum is in reality a series of closely spaced lines which extend 

to the near infrared up to about 1.3 microns. 

A number of investigators have examined these structures using both 

absorption and emission studies. The absorption of N0 9 at low temperatures 

"l A 

(liquid helium) was studied by Robinson 1 . The structure at this tempera- 
ture remains complex, indicating extensive band structures. The visible 

NO bands have been observed as chemiluminescence in the reaction between 
2 

NO and atomic oxygen and it is probable that the air afterglow is a 

pseudocontinuum produced by the many overlapping bands in this reaction. 

For NO^ produced by exciting lines photochemically, the lifetime of the 

excited state is found to be about 44 microseconds. As noted by Hertzberg 

this lifetime is two orders of magnitude longer than the value calculated 

20,21 

using integrated absorption coefficient data. Douglas suggests that this 

difference is connected with strong perturbations of the excited state by 

the numerous vibrational levels of the ground state. More recently, 

22 

Douglas and Huber have analyzed a progression of bands among the visible 

absorption bands with a spacing of 880 cm''*'. These bands consist of simple 

P, Q, and R branches and are assigned to the K' = 0 *— ►K' 1 = 1 sub-bands of 
2 2 

a — ► A^ electronic transition. 

The assignment of the experimental spectra obtained in this program 
to a fluorescence of NOg was made by a comparison of the spectral flates 
with absorption and fluorescent emission studies reported in the literature. 
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Such a comparison shows that many of the bands present are the same 
as those reported by Douglas and by Douglas and Huber* In the region 
between 4100 and 480 q£ a progress ion* involving an 880 cm ^ interval 
is seen. Band heads at 4020 and at 4260 correspond to those re- 
ported above . 

A difficulty in these assignments in the relatively sharp 
features of the spectra. Instead of a relatively diffuse continuum 
a series of rather well defined lines is recorded. One explanation 
for this effect may be the production of an excited metastable NOg 

pQ 

state directly in. a one-step reaction. As reported by Broida 
the spectra obtained using an emission produced when atomic oxygen 
and nitric oxide are mixed has a definite structure, and therefore 
direct reaction 

WO + 0— * N0 2 + hp 

In order to explain these results it was necessary to develop a 
kinetic reaction sequence involving the following sequence of re- 
actions in which the energized NO^ formed via the reaction of NO 
and 0 is collisionally stabilized to form a metastable N0 2 * molecule 
in an excited electronic state. The NOg* can then either radiate or 
be collisionally stabilized to form NOg in the ground state. The 
fluorescence then is exhibited as a discrete spectrum in the bands 
observed. This process can be shown by reference to a potential 
energy diagram taken from Broida and shown in Figure 42 . 

The reaction between NO ( 2 n ) and 0 (~T) results in a 
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repulsive state B or a bound state C. Radiation from C to A directly 
would be expected to result in a relatively ill-defined spectrum be- 
cause of the shallow potential minimum exhibited by the state C. 
However, the transition of C to B can occur by a resonance transfer 
process . Molecules formed in B can radiate to the ground state A and 
because the potential minimum is deep, a complex structure can result. 
The spectra thus observed corresponds to the experimental spectra 
obtained in this program. The major difference is an increased defi- 
nition of major structure. 

Under the present experimental tests it is unlikely that a 
direct reaction between NO and 0 occurs. There is no evidence for 
the presence of appreciable quantities of atomic oxygen in the spectra, 
although further work' should be performed in the ultraviolet regions 
to verify this point. Under the present reaction conditions a direct 
formation of the upper electronic state of NOg (B) in a bimolecualr 
collision between NO and NO^ is much more likely. In this case, 
reactions (l) and (2) are combined in one step, resulting in an oxy- 
gen transfer reaction. The subsequent reaction sequence remains 
unchanged and involves either a radiation or collisional deactivation. 
Because the NO 2 (b) is formed directly, the reaction NO 2 (b ) — *■ NO 2 (c) 
must be relatively slow compared to the probability of radiation to the 
ground state, thus accounting for the clarity of the spectra. 

An estimate of the spectroscopic value for the heat of reaction 
of NOg — ► NO + 0 (all ground electronic states) can be obtained from 
the observed short wavelength cutoff at 4060 +2QA.. Examination of 
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the weaker emission bands indicate the presence of a number of indi- 
vidual band or atomic lines . In addition to the presence of Na and 
Li in the spectra, other components which have been tentatively 
identified are CH, OH, the further positive bands of Ng, CN, NO, 
weak NH 2 , and others. These species would all be expected to result 
in gas phase reactions at high temperatures, and therefore their 
presence is predicted. 

Under most circumstances these bands would be much greater in 
intensity and thus their relatively weak intensity must be explained. 
Because of the relative weak intensities observed, it is quite likely 
that these molecular species are present in relatively low concen- 
trations. Based on calculated integrated emission probability func- 
tions the resultant spectra should be much stronger. The most pro- 
bable explanation for the lack of intensity is that the concentration 
levels are low, corresponding to a nonequilibrium distribution of 
propellants. Therefore, it is inferred that kinetic processes in- 
volving the vaporization of the fuel component may be important. 

This inference is supported by other lines of evidence including 
direct photographic documentation of the startup of the chemical 
reactor. 

Perhaps as important to the above argument as the species ob- 
served is the absence of major expected emission bands. There is 
little evidence that the NHg alpha emission bands are present or that 

hydroxyl radical is present in the reaction mixture. Although the 

\ 

strongest emission band for hydroxyl emission is the 3064.R band which 
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is interfered with by the fore optics of the spectrograph, in the 
tests using a Cini quartz prism spectrograph no OH was found. This 
also tends to confirm the possibility that there is a high degree of 
nonequilibrim. 

5 . 4.5 HNO in Emission 

A typical "real time" emission spectrum from MMH/NTO engine 
firing is presented in Figure 43 . Analysis of this spectrum indi- 
cates that the HNO (100 - 000 ) transition at 6 l 7 Clft may be present in 

emission. The assignment of this band is based on the work of Cashion 
24 

and Polanyi . Unfortunately, because of a dropoff in the film 

o 

recording sensitivity above the 65OQA. region, confirmation of the 
(0,1,0) *—♦ (0,0,0)' 1 and (0,0,0) '— *■ (0,0,0) " bands at 6930ft and * 
768($ respectively, could not be made. Additional data is necessary 
in the near infrared region in order to remove possible ambiguities 
in the band assignment. The presence of this band correlates with 
infrared absorption data obtained from engine residues and may indi- 
cate that metastable intermediate products condense on the engine 
walls. The lifetime of HNO in the excited state is long relative 
to the pulse firing duration and has been estimated to be 0.1 second 
by several investigators . This is much longer than the total firing 
duration of our present engine. 

5 - 4.6 Physical Characterization of the Residues 
. The liquid residues formed in a pulse mode engine had a con- 
sistency corresponding to a medium to high viscosity liquid. The 
material was yellowish-brown, which after a period of several days 



89 



4 r/t/E 


6300 6200 6100 6000 5900 5800 5700 



Figure 43. Hydrazine -Nitrogen Tetroxide Spectra 
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turned into a dark aquamarine. The specific gravity was 1.2 for a 
residue formed from an MMH/NTO test. The residues under vacuum had 
a frothy characteristic and were observed to outgas or bubble at a 
relatively fast rate. A photomicrograph of a typical residue is 
shown in Figure 44 . As can be observed, a large number of gas bub- 
bles is present. The largest bubble diameter (at ambient pressure) 
was 1.2 mm. A plot of the distribution of sizes in a typical area 
is given in Figure 45 . It can be seen from this plot there is a 
relatively large number of smaller bubbles in the mixture, but fre- 
quency of occurrence is relatively constant in the range 0.2 to 0.5 
mm. The fact that the bubbles coalesce into each other accounts for 
the smaller number at the larger diameter. The maximum bubble di-' 
ameter is limited by the depth of the film thickness . 

5.4.7 Pressurization Studies for N 2 °4 

A series of time/pressurization curves were developed for NTO 
injected into air at various start pressures . The experimental 
results are shown in Figures 46 through 47 . As expected the rate of 
pressurization is dependent on the initial start pressure and is much 
more rapid at low initial pressures, decreasing with increasing am- 
bient pressure. A comparison of the rate of pressurization for a 

2 

NpO^ pressurization calculated by the CSC program was made and a 
typical run presented below. 
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Figure 44. Photomicrograph of Propellant Residue 
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NUMBER OF BUBBLES 







Fig. 46 Pressurization Rate Histories for a N 2 0^ (l) Jet Stream Under Simulated 
Space Start Condition. 

(Top trace: Pressure, Bottom trace: Flow, Time base: 5 ms/cm, N p 0, flow 
volume: 1.5 cc, Transducer pressure: 50 lb/cm) ^ 


94 


Run 26. Start pressure: 9«0 psia Run 27- Start pressure: 12.0 psia 










Fig. 47 Pressurization Rate Histories for a N 0^ (l) Jet Stream Under Simulated 
Space Start Condition. 

(Top trace: Pressure, Bottom trace: Flow, Time base: 5 ms/cm, NO, flow 
volume: 1.5 cc, Transducer pressure: 50 lb/cm) 2 ^ 
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Run 21. Start pressure: 20 mm Hg. Run 22. Start pressure: 20 mm Hg 






Pressure 


Time 

Experimental 

Calculated 

0 ms 

0 psia 

0 psia 

1 ms 

10 psia 

3 psia 

2 ms 

17 ps ia 

10 psia 


As can be readily observed there is a difference of a factor of two or 
three between the rate of pressurization experimentally determined and 
the rate as calculated in the preignition time regime using the Seamans 
program. This study was conducted in order partly to provide a better 
understanding of the vaporization rates as they affect the ignition and 
partly to document the phenomena using the unique flow time analytical 
capability of the rocket engine developed on this program. 

5.5 Ammonia/Hydrazine Mixtures 

One of the obvious problems was that the one-on-one injector used 
in this study did not produce the same droplet size distribution for the 
hydrazine fuels and the nitrogen tetroxide oxidizer. This fact makes 
any interpretation using analytical models involving uniform droplet 
sizes (e.g., Seamans 2 program) extremely questionable. Two solutions to 
this problem of increasing the breakup of fuel sprays were investigated. 
One solution involved the use of a chemical additive and the second in- 
volved the use of heated fuels in which the fuels were heated up to a 
temperature at which the vapor pressure of the hot hydrazine was equal 
to the vapor pressure of nitrogen tetroxide at normal temperature 
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injection conditions. Because the use of chemical additives imply a 
reduction of performance it was desirable to illustrate the principle 
by using an "additive" which would in itself be a high performance 
propellant fuel. Ammonia (NH-^was chosen for this purpose and mix- 
tures of ammonia/hydrazine were injected into the space chamber using 
the injector block assembly previously described. Typical results of 
this test are shown in Figures 48 a,b,c. Figure 48a shows ^0^ in-' 
jected at about 15 mm Hg, Figure 48b shows injected at about 

10 mm Hg,and Figure 48c shows a 52/48 percent by weight mixture of 
MH^lJ/N H^(l) injected at 10 mm Hg. As can be readily observed, 
there is little breakup of the N^H^, a very rapid breakup of N^O^ 
with a characteristic fan at the injector face and a reasonably good 
dispersion of the NH^/N^H^ propellant under these conditions. It was 
concluded that if it became necessary to increase the dispersion of 
liquid hydrazine fuels it could be accomplished by changing to a 
'NH-/N 2 H 4 mixture. It is also noted that theoretically it should be 
possible to make substantial changes in the rate of propellant stream 
breakup by modifying the surface tension and that the use of surfactants 
(surface active agents) used in the stabilization of gels, etc., should 
be investigated. There is a problem of the transport of the surfactant 
molecules to the droplet surface under the very short times associated 
with the injection process. 

5.6 Heated Propellant Dispersion Studies 

The above argument is based on the concept of changing the vapor 
pressure of the fuel to enable a sure rapid volatilization and breakup 
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a * N 2 0 4 at 1 5mm Hg 




Figure 48. 


Comparison of Spray Patterns for N,0 
52/48 Mixture of N H 3 (i)/N 2 H 4 (f) 
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of the jet stream. Theoretically, it should he possible to accomplish 
the same effect by heating the propellant fuel to a temperature corre- 
sponding to the temperature at which the vapor pressure of the fuel jet 
is approximately equal to the vapor pressure of the oxidizer at ambient 
temperature. A series of hot fuel tests were conducted to verify this 
hypothesis. The results of this series of tests are shown in Figure 49 . 
As can be readily observed there is only a very small change in the 
dispersion of the jet of hydrazine at 165 F. Therefore, although it 
may be beneficial to heat the propellant in order, for example, to 
increase the reaction rates, the heating apparently does not make a 
substantial change in the rate of jet disintegration under these ex- 
perimental circumstances. 

5.7 Real Time Pressure and Visible Emission Tests for a 
NgH^/N^O^ Propellant 

5.7.I Pressure Studies 

Experimental pressure data were obtained for a rocket 

test firing. The test conditions included startup using a clean engine 
at a simulated altitude greater than 90,000 feet. During this test one 
of the two l/4 in. quartz viewing ports was blown out. The fracture 
pattern indicated destruction probably occurred via a high overpressure 
rather than a detonation at the window site. Propellant flow was 
characterized by a 2 msec fuel lead and total flow time was 

approximately 20 msec. The pressure-time trace indicated a series 
of three pressure spikes each in excess of 700 psi. The individual 
pressure spikes are separated by 3 to 4 msec indicating a low 
frequency chugging at 250-300 cps . There was a one-to-one correspondence 
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a. N 2 H 4 Dispersion at 125°F 



b. N 2 H 4 Dispersion at 165°F 


Figure 49. 


Effect of Fuel Temperature on Hydrazine Jet Spray Pattern 


100 




between these low frequency pressure fluctuations and observed vari- 
ations in the fuel and oxidizer flow rates . 

A composite pressure-time history of the test is given in Figure 50 
through 52. This history covers the time immediately prior to ignition 
to the onset of the post combustion time period. In this case the 
time scale was expanded to 100 jusec - cm - "' - in order to be able to 
correlate pressure fluctuations with emission spectra obtained from the 
W/S scanning spectrophotometer operated at a 1.25 msec sweep rate. 

Immediately after ignition a base line shift occurs for the 
pressure transducer output. This shift requires a correction of 300 
psi to be applied to the data. It was believed that this shift was 
caused by the thermal shock of the piezoelectric crystal. In order, to 
verify this hypothesis a series of tests were performed in which radiant 
energy emanating from an A-l flashbulb was directed onto the surface. 

A direct correlation was made between the flash duration and the base 
line deflection of the transducer under conditions involving no change 
in pressure. Although the rise time of the flashbulb is slower than 
the observed rise time of the spectral radiation and pressure in the 
2-dimensional rocket engine, the grey body emissivity of the bulb is 
about the same as that expected from the rocket combustion. From this 
evidence it was concluded that the base line shift is caused by the 
influx of radiant energy onto the surface of the transducer and hot 
to other possible sources of instrumental error such as point-contact 
inductances . 

Although the physical and chemical processes occuring in the 
thrust chamber are very complex, the form of the correction curve to 
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be applied to the transducer output can be obtained by observing that 
immediately following the high pressure transient a series of rare- 
faction waves occur. Evidence for this effect is indicated by the 
voltage output falling below the original base line. This envelope 
of points provide a minimum value to the rate of shift of the base 
line caused by thermal shock, since it can be shown that the value of 
the pressure during this rarefaction cannot be less than zero. Using 
this envelope of points the actual chamber pressure can be estimated. 

The gross pressure-time history indicates a series of three high 
pressure transients. It is probable that each of these three transients 
represent a separate ignition event since direct propellant flow mea- 
surements indicate that the chamber pressures are high enough to 
reverse the flow in the injector manifold and propellant tanks. The 
result of flow reversal caused by a preceding pressure spike is that 
the next overpressure spike results from propellant flow conditions 
with zero fuel lead or lag. Therefore a number of similarities should 

exist in the detailed pressure-time history of each overpressure spike. 

I 

The structure of each "ignition spike" appears to be similar. 

A gross description indicates that there are four regions of interest. 
The first region consists of the pre ignition time period in which no 
increase in pressure is observed (within experimental error). The 
second region is the initial buildup of pressure in the engine. This 
occurs in time period of 30-50 H sec during which time the chamber 
pressure rises from a very low value to a pressure level between 75 
and 150 psi. The shape of the pressure curve is expotential in form 
indicating the onset of a "runaway" temperature excursion. This is 
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to... 


typical of combustion processes accelerating from a normal deflagration 
to a detonation. The third region occurs over a time' period of about 
20 j* sec. During this time period a closely spaced series of high 

pressure transients are observed. The pressure levels reached during 
this time period are in excess of 1000 psi with rise times on the order 
of between 1 to 2 (i sec Because the response of the 603^ trans- 

ducer is limited to 1 ^ see it is possible that the actual pressure 
rise times may be considerably shorter. The fourth time region 
probably involves a "more or less steady state" deflagration of the 
propellants which results in a pressure decay from about 700-800 psi 
to a low steady state chamber pressure . 

A preliminary analysis associated with the third region of 
interest indicates that there are at least three superimposed wave 

5 4 4 

forms with individual frequencies of 4x10 , 4x10 and 2x10 cps 
respectively. The high frequency component corresponds to the crystal 
ring frequency of the 603A transducer. The low frequency component 
can be interpreted as a standing wave in either the longitudinal or 
lateral resonance modes of the thrust chamber. The component under- 
goes a rapid attenuation. If this assignment , is correct then an 
average chamber temperature in the neighborhood of 3000°K can be 
derived implicitly from sonic velocity considerations. Two alternate 
assignments can be made to account for the intermediate frequency wave 
form. The first alternative involves the excitation of a standing 
wave oscillating in the vertical mode (short path length) of the 
thrust chamber. This assignment then implies a localized temperature 
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between 5000-6000°K. A second alternative is to assign the intermediate 
oscillation to the excitation of a standing wave in the thrust chamber 
walls. The implied sonic velocity of 20,000 fps is similar to that 
expected for an aluminum thrust chamber. This latter assignment re- 
quires a relatively close coupling of the thrust chamber with the trans- 
verse oscillation mode of the piezioelectric crystal. Experiments in 
which an explosive actuation of the normal resonance frequencies of the 
thrust chamber are planned for the future in order to aid in the inter- 
pretation of the transducer output signals in this intermediate range. 

5.7-2. Spectral Emission Studies in the Visible Region 

5. 7. 2.1 Analytical - The Effect of Pressure on High Temperature 
Gaseous Emissivities 

As described previously, strong pressure fluctuations are observed 
during the ignition and combustion of the pulse mode rocket 

engines. In order to interpret the spectral emission obtained during 
engine test firing it is necessary to understand the relation between 
molecular emissivities and pressure. A preliminary analysis of the 
dependence of emissivity on pressure is outlined in this section. The 
treatment presented herein closely follows that of Penner. 

The effective emissivity ^ of an individual species in a high 


temperature gas is given by the expression 

e 


4.76 x 10 7 /hc\ 5 4 

ef f , i 26.4 YkTj % SXp 


^~ hc 


( 1 ) 


where 


h = Plancks constant 

k = Boltzmann constant 

T = absolute temperature 
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the velocity of light 
the fundamental vibrational frequency 
path length 

the absorption oscillator strength 

the ratio of the number of species i molecules to 
the total number contained in the system 

the density of the system relative to its density 
at STP conditions. 

The spectral distribution of radiant energy for each species is 

given by 

d6 eff,i = 4q) 2 diju (2) 

du) eff,i uo 4 

o 

For each molecular species (e.g., NO) it is necessary to sum over each band 
that can contribute to the radiant energy. 

For NO, the G and y bands are of major importance in the visible region 
These bands correspond to the 

B 2 T X 2 7T and A 2 E + -**X 2 Jf ( 3 ) 

transitions respectively and 


UU 


o 

6 

f 

♦ 


P/P, 


e. . = E £ e . (4) 

tot j i eff,i v ' 

From this sum it is possible to determine both the shape of the 


emissivity curve as well as the absolute magnitude of the total energy flux 
in any specific interval of wavelength. 

Pressure effects on emissivity occur from interactions of’ the last 
four parameters of equation(l). As a first approximation it can be assumed 
that the major effect will result from variations in p /p . This assump- 
tion results in a linear dependence of effective emissivity on pressure. 
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The effect of ignoring the f parameter is to legislate that the pressure 
fluctuation frequencies are much faster than the characteristic kinetic 
relaxation time of the system. It is probable that this parameter can be 
neglected for the intermediate frequency pressure oscillations but that 
substantial errors may exist in correlating emissivity data for the low 
frequency pressure spikes. The effects of pressure on the 6 and f 
components of emissivity can be estimated via considerations of optical 
depth and Doppler line shifts . 

5 . 7 • 2 . 2 . Direct Determination of Temperature from Experimental 

Spectral Data 

5. 7- 2. 2.1 Theory 

The two-color technique for the determination of temperatures from 
spectroscopic data requires an experimental determination of the radiant 
intensities emitted from a given flame in two well-defined wavelength 
regions. The measured radiant intensities define two brightness temper- 
atures and a color temperature. The brightness temperature T br at the 
wavelength X is defined as the temperature at which the emitter under study 
emits the same spectral intensity of radiation as a blackbody. The color 
temperature T c of an emitter is defined as the temperature at which a 
blackbody emits radiation having the same ratio of radiances at the wave- 
length X and X 2 as the fitter under study. 

R°(X 1 ,T C ) = e pU^T) R°(X I3 T) 

r°(x 2 ,t c )' £ 2 (x 2 ,t) r°(x 15 t) 

This expression assumes that the two brightness temperatures are the 
same and the emissivities of the species in the flame zone are invariant 
over two chosen wavelength regions. If the emissivities are not the same, 
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then independently determined emissivity data must be used in order to 
calculate the true flame temperature from the observed consideration of 
the individual emitting species providing the relative concentrations are 
know, but the procedures are relatively complex. 

As a first approximation the assumption of constant emissivity will 
be utilized in this treatment. The rate of emission of radiant energy 
from a blackbody, i.e., the emitted radiant flux, is then given by Planck's 
equation. For sufficiently small values of the product of wavelength and 
temperature, Wien's equation is applicable. 

R°(X,T)dX=% exo(-— )dX 
X 5 X7 

Combining equations (l) and (2) the relation between the color temperature 
and the true temperature can be derived. 

1 1 - £n(ei/e 2 ) 

T* T c C 2 [(lA 1 ) - UM 2 )3 

This relation indicates that the color temperature will be equal to the 
true temperature, independently of the choice of X^, andX 2 if the emit- 
ting source is either a blackbody (s =e 2 =l) or a greybody ( 6^ =e 2 <i). 

j 

A direct test can be made for the presence of a black or greybody by 
determining whether the observed color temperatures are independent of 
wavelength . 

5. 7-2. 2. 2 Experimental 

A direct determination of temperature using a two-color spectrographic 
technique was made for a propellant combination. Since each 

scan required 1.25 msec, the time resolution of the temperature measure- 
ment corresponded to this interval. A series of trial calculations were 
made to establish which spectral region exhibited greybody emission charac- 
teristics. Based on these trial calculations, a color temperature derived 
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from the relative radiant emissivities at 540oX and 500 oX was selected, 
based on the observation that the emissivity in this region was found to 
be relatively independent of wavelength. The derived temperatures for 
Test 12 are given in Table II. As indicated on several scans, there is 
an undefined temperature corresponding to large ignition spike over- 
pressures which occurred during the time interval. 

Based on this data it is interesting to note that temperature rises 
between 400 to 1000°K-msec ^ are observed immediately prior to an ignition 
spike. Also, it is noted that after complete propellant shutdown there 
is continued burning and that the rate of temperature decrease is roughly 
linear at about lOC^K-msec - ^. A summary of the flow conditions for the 
burn include a short fuel lead (3 msec) with a nominal O/F mixture ratio 
of 1.8. The termination of the oxidizer flow occurred during the third 
ignition spike and the fuel flow terminated approximately 1.0 msec later. 
The induction time required for the first ignition spike was approximately 
7 msec respectively. A space simulation altitude in excess of 100,000 ft 
was used for this test. 
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TABLE II 


Spectroscopic Temperature for N^H^/NgO^ 


Test 12 



1.25 
2.50 

3.75 

5.00 

6.25 
7.50 

8.75 

10.00 


11.25 

12.50 

13.85 

15.00 



3100 

3000 

2800 

2675 

2450 


100 

125 

225 


Comment 

First Igniton Spike 


Second Ignition Spike 

Third Ignition Spike 

Termination of 
Propellant Flow 
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VI. DISCUSSION 


6.1 Introduction 

Based. on an extensive hot firing program, several engineering 
correlations had been developed with regard to pressure spiking in 
pulse-mode NTO/NgH^ rocket engines. These correlations had been 
established by the following: (l) the catastrophic overpressures 

were a function of previous engine test history; (2) the probability 
of the overpressures was increased for engines which were fired in an 
upward attitude; (3) the probability of large ignition overpressures 
was increased as the ambient temperature of the rocket engine decreased; 
and (4) the overpressure severity increased as the ambient pressure 
in the thrust chamber decreased. These correlations indicated the 
existence of a causal relationship between engine residues and the 
high pressure spiking phenomena. 

Additional evidence for such relationships was available from 
high speed photographic studies of the rocket plume performed at NASA- 
Hbuston. In this study it was observed that a liquid residue was 
formed in the engine during actual test firing and that residue was 
observed to drop from the lip of the exhaust nozzle. This observation 
was unexpected since thermodynamic calculations for the engine indicated 
a chamber temperature well above the boiling point for any of the con- 
sidered products of reaction. The fact that this material was a liquid 
under these extreme environmental conditions indicated that the material 
was neither the raw propellant nor products of reaction such as water 
produced in the process of combustion. 
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Preliminary tests at Marquardt indicated that this material 
might consist of a combination of fuel and fuel nitrate. These 

p 15 

chemical tests at Marquardt were performed by showing that 

a liquid phase reaction between and HNO^ resulted in an 

"energy rich" product which underwent an explosion when confined 
in a sealed capsule. A direct chemical analysis of the residue 
for nitrate ion was positive. 

In a subsequent research program at the Bureau of Mines, 

q U s 

Perlee ’ ’demonstrated that the compounds were, in fact, hydra- 
zinium nitrates. This demonstration was made by synthesizing 
hydrazinium nitrate, unsymmetrical dimethyl hydrazinium nitrate, 
and monomethyl hydrazinium nitrate by means of standard chemical 
synthesis, and then comparing the infrared spectra of these 
compounds with the infrared spectra of actual engine residues. 

1 2 

This analysis was correlated by results obtained by Seamans ’ 
who also performed an elemental chemical composition analysis 
that indicated a mixture of hydrazinium nitrate admixed with raw 

j 

fuel. Based on this evidence it was generally agreed: (l) that 
hydrazinium nitrate was formed; and (2) that it acted as a primary 
explosive and was responsible for the detonation. 

In addition to the hydrazinium nitrate, Dauerman u reported 
the presence of a short lived azide peak at 4.9 f* • This transient 
peak occurred in the condensate formed from the direct reaction 
between nitrogen dioxide and hydrazine under closely controlled 
laboratory conditions. Support for this possibility was also found 
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from data obtained by B. Clark* at Bell Aerosystems and by H. 0. 
Weiss^ at Dynamic Science. In this work it was necessary to 
freeze the residues at LN2 temperature in order to obtain the k.9yU 
absorption peak. Infrared studies at 4.^tt conducted by Perlee, 

Seamans, , and at Edwards Air Force Base (unreported)' were negative 
for azide. Therefore, there was a difference in experimental results 
with regard to the existence of an azide in the propellant residues. 

A major reason for being interested in identifying azide deriv- 
atives as components in propellant residues was that studies by Perlee 
indicated that shock overpressures greater than 22 kbar were 
required for the initiation of hydrazinium nitrate. Such pressures 
had not been observed during hot firings and therefore it was believed 
that a shock sensitive chemical initiator had to be present in the 
residues. Because of their inherent instability and their relative 
ease of chemical formation the most obvious candidates would be the 
azides . 

A difficulty with the identification of azide intermediates was 
the nature of the azide peak, which as observed by Dauerman had a 
"lifetime" of about 8.5 sec. In a pulse-mode operation a typical 

firing cycle will involve a series of short bursts 13 msec on, 87 msec off. 
The lifetime of the azide peak then was long relative to the "off time." 
However, in a number of instances high overpressures occurred upon 
engine starts which were initiated long after a previous test. The 
longest documented occurrence of this type was about 1.5 hours after a 
previous test. Partly because of this test and partly because of the 

*Personal Communication. 
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lack of agreement between several investigators as to the presence of 
hydrogen azide in the residue material, an alternate mechanism was 
invoked by Perlee to explain the initiation. This alternative was the 
hypothesis that the initiation was caused by a rapid cavitational 
collapse of a void in the residue created by a shock reflection at the 
chamber wall. The shock was assumed to have been generated during a 
relatively "low pressure" ignition overpressure. This possibility has 
not been disproven as yet and may, in fact, be a source of initiation 
for certain engine geometrical configurations. 

Seamans- 1 -^ observed another characteristic effect associated 
with the residue. In a series of Trauzl block tests he observed that 
about fifty percent of the time the lead test block was completely 
destroyed. A similar occurrence took place during hot firing programs 
at NASA-Hbuston in a test engine. This high brisance of the residue 
was difficult to explain for two reasons. First, the residue consisted 
principally of a mixture of water, hydrazine, and hydrazinium nitrate. 
Thermochemically this mixture was a very poor explosive which would be 
expected to have a brisance very much less than nitroglycerin which was 
used as a comparative standard. Secondly, the more sensitive materials 
in these residues, which could conceivably cause this high brisance, 
were present only for relatively short periods of time. Also, there 
was not direct evidence that these intermediate compounds were present 
in the residues analyzed by Seamans . 

Based on the above survey, three major questions arose. The 
first question was concerned with the chemical composition of the 
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residues, i.e., to what extent were chemical intermediates formed 
and what was their chemical characterization. The second question 
was concerned with the extent of chemical disequilibrium in the 
preignition and combustion phases of operation. Generally it was 
assumed that pseudoequilibrium conditions prevailed throughout the 
actual engine firing cycle and that the bulk of the residues were 
formed by the condensation of raw propellants on the chamber walls 
after engine shutdown. The source of these propellants was the 
residual material which remained in either the injector manifold or 
the fill lines after the engine valves were closed. The third 
question which had to be answered was how to explain the high bri- 
sance characteristic of the detonation and a corollary question of 
why the phenomena had occurred only during high altitude simulation 
tests and had not been observed for the sea level testing incurred 
during previous engine development programs. The following discussion 
outlines test results obtained during the course of this research 
program and suggests a physical mechanism based on these results which 
can explain all of the major phenomena observed. 

6.2 Chemical Characterization of the Engine Residues 

The chemical analysis of propellant residues has been described 
in sections IV and V. Infrared absorption studies of these 

t * • 

residues confirm the results of Perlee^’^’^, These studies indicate 
that the primary composition of the residue consists of a mixture of 
the hydrazine and hydrazine nitrate. A direct correlation can be 
made with the infrared spectra of the hydrazinium nitrate derivatives 
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synthesized by Perlee and the infrared spectra obtained from actual 
residue material. There is little doubt that the residues formed 
consist primarily of these compounds together with combustion pro- 
ducts (e.g. water). 

There is also direct infrared evidence for the presence of 
hydrogen azide, which tends to support the experimental results of 
Dauerman. The spectrum obtained in this program shows a much weaker 
absorption pattern, indicating only a small quantity of material. 

This is in contrast to the strong absorption peak observed by Dauerman. 
The most probable reason for this discrepancy is that under the ex- 
perimental conditions used by Dauerman the early liquid - liquid phase 
reaction was observed and the volatility of HK^ was sufficiently high 
to result in the loss of a large percentage of the material after a 
substantial period of time. 

The presence of HN3 and/or other azides was also confirmed by 
mass spectral analysis. For the mass spectral analysis the sampling 
procedure involved an "in vacuo" degasification of the residue. 

Samples of the effluent gas produced by bubble formation were then 
obtained by allowing the sample to froth for several hours . This pro- 
cedure thus allows an analysis of the gases contained in the bubbles 
produced by the slow liquid - liquid reaction. By this method these 
bubbles were found to contain HN3, triazene, and methyl triazene for 
residues obtained from an MMH/NTO propellant combination. Other com- 
ponents such as hydrogen and nitrogen were probably also present. 
Quantitative analyses of the relative mole fraction of individual 
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species were not possible without an extensive study of individual 
breakdown patterns in the mass spectrograph ionization chamber. 

Except for the hydraz inium nitrates, the compounds were present in 
very small quantities. Assuming that these materials are distri- 
buted more or less evenly throughout the liquid it is doubtful that 
they would participate as a liquid phase initiator for the primary 
hydra.zinium nitrate explosive. 

From these results we conclude that the residues contain a 
relatively large number of explosion- sensitive compounds. However, 
these compounds are present only in trace quantities. 

6.3 Mechanism and Theory of Explosive Initiation in Rocket 
Engine Residues 

6.3.1 Mechanism 

The residue formed from an NTO/N 2 H 4 propellant consists of 
a mixture of fuel and fuel nitrate together with small amounts of 
detonation- sensitive compounds. These residues are also character- 
ized as high brisance explosives, because extensive damage results 
when they are initiated. In rocket tests this initiation has occurred 
only during high altitude simulation testing. It is necessary, there- 
fore, to explain why a thermochemically unbalanced and relatively poor 
explosive can exhibit high brisance detonation behavior, and why this 
behavior occurs only during the simulated high altitude testing,. 

The theoretical approach used to describe these processes 
assumes that the engine pressurization processes cause an adiabatic 
compression of small bubbles that are formed in the liquid residues 
by slow "liquid phase" reactions. The pressurization process then 
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results in the development of a localized "hot spot." This "hot 
spot" may be inert or chemically reactive. If it is chemically 
reactive then its temperature involves a balance between the rate 
of heat dissipation caused by convective and conductive heat fluxes. 
The critical condition for explosion occurs when the self- heating 
within the spot exceeds the rate at which heat can be dissipated. 

Based on the experimental evidence obtained during the course 
of this program, it is assumed that the bubbles consist either of 
pure hydrogen azide or hydrogen azide diluted by other components. 
Since the concentrations of these components is not known it is 
assumed that the dilutent may be described as an inert substance, 
e.g., nitrogen. 

6.3.2 Hot Spot Theory 

A number of theoretical formulations have been developed to 
describe the critical conditions for the occurrence of an exponential 
temperature rise resulting from a chemically reactive hot spot con- 
tained in an explosive material. These formulations are 
based on a solution to the general differential equation. 


f>c|i = k/^-,+ 


+ 


Qfe' 


e/rt 


where K, P , and C are the thermal conductivity, density, and specific 
heat respectively and are assumed to be temperature independent.. P is 
a steric factor which assumes the values 0, 1, or 2 for a slab, cylin- 
der or sphere, and Q, E, and f are the heat released per unit volume, 
the activation energy, and the Arrhenius pre-exponential factor, 
respectively. 
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The initial "boundary conditions are assumed to be : 


T = Ti x < A 

T s To x > A 

where A is the half width or radius of the hot spot. 

With these boundary conditions, the dimensionless variables 
0 and 8 can be formulated: 

e z [e/rt’] [ T 1- T .] (6.3-2) 


8 = 


Q f E A 2 
RT?K 


-E/RT 


(6.3-3) 


This mathematical formulation neglects loss of reactant, change of 
thermal properties with temperature, differences in the thermal 
properties of the hot spot and the material outside of it, changes 
in phase caused by thermal heating and reactions other than single 
step reactions controlled by temperature. For a given 9 these vari 
ables assume critical values 0 C and 8 t as a condition of explosion. 

The relations for the critical condition have been examined 
by previous investigators and a summary of the correlations between 
§ e and 0 C for hot spot spheres are listed in Table III. From these 
correlations the critical radius can be evaluated at any specified 
temperature by 

A = [ RT, 2 k 1 [of E e" E/RT i 1 [*»] 

L JL (6.3-4) 


where f(9 0 ) is defined from Table III values. 
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TABLE HI 


CRITICAL COMPARISON OF CRITERIA 
FOR SPHERICAL HOT SPOTS 


Boddington 


Friedman 


Zinn 


Merzhanov 


Rideal and Robinson 


Thomas 



6 C = 4. 7 6 0 
6 C = 25 

0 . 6 

6 C = 12. 1 In 9 0 

6 c = 2 0 q 

r 25 

§ r = — - T 

I0tt2 50tt 

+ e 0 (1 + P) + 0 O 2 (1 + P) 2 
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The critical induction time, r c , can be evaluated using 


t .p± A ( 6 . 3 - 5 ) 

Sc* 

where A c corresponds to the critical radius and c, K, and P are 
evaluated at the hot spot temperature, T p , attained after an adia- 
batic compression. The thermal conductivity K was calculated using 
the Eucken approximation for polyatomic molecules 



( 6 . 3 - 6 ) 


In the viscosity equation a represents a Lennard- Jones character- 
istic cross section and A is a temperature- dependent Stockmayer 
potential function. C p is calculated using the sum of the trans- 
lational, rotational, and vibrational partition functions and the 
viscosity m is approximated using 


,^-s/mT 

M = 2,669 xlO — 


( 6 . 3 - 7 ) 


A M 


In order to apply these theories it is necessary to estimate 


the values to be assigned to the activation energy E and pre-exponen- 
tial factor f. For a unimolecular hydrazoic a.cid decomposition a 
value of 10' J was assigned to the pre-exponential factor. The acti- 
vation energy was estimated at a value of 40 kcals/mole. This value 
for activation energy was based on the enthalpy change for the reaction 


HN (V)-* NH ( X A ) + N q ( X 2 + ) 


(6.3-8a) 


• o \ ** / " ^ — / " 2 

In this reaction it is noted that dissociation into the NH ( 'A ) 
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state is assumed rather than dissociation via either of two other 
possible routes: 

hn 3 ( 1 A’ ) — *NH ( 3 s") + N 2 ( ls g + ) (6.3-8b) 

hn^ (V) — *n 3 ( 2 n g ) + H ( 2 S g ) (6.3-8c) 

Reaction (6.3-8c) requires an activation energy of at least 85 
kcals/mole and is therefore ruled out. Reaction (6.3-8b) although 
energetically most favorable, violates the Wigner spin correlation 
rules and is quantum mechanically "forbidden." 

The energy of the *A potential curve above that for £ was 
taken as 1.3 electron volts although it should be noted that uncer- 
tainties in the correlation energy would imply only a value less 
than 2.0 ev . Values for the enthalpy of NH were taken from the JANA.F 
Thermo chemical Tables and the heat of formation of HN^ was taken to 
be +73 kcals/mole. Using the above values as minimum, actuation 
energy of 40 kcals/mole was estimated for reaction (6.3-8a). Together 
with the correlations outlined in Table II for the several hot spot 
theories, values of d c , A c and t c were calculated. 

A graphical analysis using these theories is given in Figures 
53 - 57 for pure HN^/l^ compositions . Of particular interest tq this 
study is the variation of critical radius and critical induction time 
with initial temperature. The initial radius corresponds to the radius 
of the collapsed bubble following an adiabatic compression. The initial 
compression process results in the heating of the bubble to an elevated 
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Figure 53 Comparison of Theoretical Predictions for Hydrogen Azide Critical Radius 
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Figure 54. Temperature Increase Caused by the Adiabatic Combus 
tion of Hydrogen Azide/Nitrogen Mixtures 
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Figure -55. Critical Initiation Times for Pure Hydrazoic Acid Gaseous Hot Spots 
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Figure 56. Adiabatic Compression Temperature Vs. Critical Postcompression Radius 
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Figure 57. Adiabatic Compression Temperature Vs. Pressure Ratio 





temperature T2, calculated using the expression 


7(T)- 1 



in which 7 (t) is the temperature- dependent ratio of specific heats. 

The value of was calculated using statistical mechnical. partition 
functions in which the translational and rotational degrees of freedom 
were considered fully active and the vibrational partition function 
was calculated using the Einstein functions and spectroscopic funda- 
mental frequencies. For HN^ the frequencies were taken as 3336, 

2140, 1275, 1150, 522, and 672 cm -1 . 

In this model the residue material during a vacuum ignition 
event is subjected to a pressure pulse corresponding to either a 
steady state combustion pressure or an ignition transient overpressure. 
For the present study the steady state combustion pressure is about 
100 psia although transducer data overpressures up to 1500 psia have 
been recorded experimentally. Assuming that steady state thrust chamber 
conditions are achieved, the pressure ratios causing an adiabatic temper- 
ature increase will vary between 50:1 and 4000:1. This corresponds to 
vehicle operations between 50,000 and 150,000 feet altitude. For pure 
HN3 the resultant compression temperatures were calculated to lie 
between 750 and 950° K. 

Using the theories previously developed for hot spot initiation 
in homogeneous explosives, calculations were made of the critical 
radius and critical time required for pure hydrazoic acid. These cal- 
culations used temperature as a variable and the results are listed in 
Figure 56 for operating pressures of 1000, 100, 6.0 and 0.3 psia 
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respectively. For operations at 50,000 feet in which ignition transients 
of 1,000 psia occur, the critical induction time is about 100 milliseconds 
and the critical postcompression radius is about 0.1 mm. For this con- 
dition the required precompression radius is estimated to be about 1.2 mu 
Since the liquid residues in the thrust chamber are between 0 - 1 mm in 
thickness we conclude that it would not be possible to initiate the ex- 
plosion at 50,000 feet. Also, since the gaseous components are complex 
mixtures and not pure HN3, it is doubtful whether heat release per unit 
volume is sufficiently high to allow initiation. This conclusion is 
supported by all the rocket test data available at the present time. 

For operations at 150,000 feet the pressure ratio is about 5000:1. 
Corresponding to this pressure ratio, the adiabatic temperature rise is 
about 950° K. Again assuming steady state conditions at the operating 
chamber pressure of 100 psia the critical time is 0.6 microseconds and 
the critical radius is 8 x 10"^ mm with an associated precompression 
critical bubble diameter of 0.2 ram. Since bubble diameters of this srize 
are observed throughout the residue it is much more probable that 
initiation will take place. It is concluded that at 150,000 feet the 
probability of initiation is high for residues having thicknesses greater 
than 0.5 mm. 

The effect of high pressure ignition transients can be estimated 
by assuming that near-equilibrium conditions are achieved during the 
initial adiabatic compression. Transducer data taken at Marquardt?^ 
indicate that pressure transients well over 2,000 psia may be reached 
during ignition. These spikes are about 50 microseconds in duration 
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and have extensive fine structure with high pressure regions on the 
order of 5 microseconds. For operations above 50,000 feet, pressure 
ratios having a minimum value of 5000:1 are reached with compression 
temperatures estimated above 900° K. Since these are essentially the 
same conditions as those found at 150,000 feet, it would be predicted 
that pressure spiking could induce initiation either by a hot spot 
mechanism during "steady state" ignition at high altitudes or during 
engine spiking at higher ambient pressure. 

The temperature rise within the reactive hot spot after ini- 
aticn can be calculated by assuming an adiabatic temperature increase 
caused by the chemical decomposition. A series of thermochemical 
equilibrium calculations has been carried out for pure hydrogen azide 
and mixtures of hydrogen azide and nitrogen. The results of these 
calculations at l4.8 psia are presented in Figure 54 for temperature 
increase and relative mole fraction. A precompression temperature of 
300° K is taken as an initial temperature for this calculation. As 
can be seen by inspection, temperatures between 2,000 and 4,000° K 
are reached following decomposition of HN q - ^ mixtures. In practice 
these temperature values would be somewhat higher due to the enthalpy 
increase caused by the initial compression. Because dissociation will 
limit the temperature increase, a direct estimation of total tempera- 
ture cannot be performed. Therefore it is assumed that these tempera- 
tures are sufficiently high to cause initiation of the primary explosive 
contained in the engine residues and that the major time delays are 
associated with the development of the reaction front in the compressed 
gaseous hot spot . 
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Example Case 


For an engine designed to operate at 100 psia chamber pressure , 
determine the critical time and critical precompression bubble radius for 
start operations of 60, 80, 100, and 120 thousand feet. For this case, 
using this theory, the following values are obtained. 


Start Altitude 

60 

80 

100 

120 

Critical Postcompression Radius (mm) 

0.5 

0.12 

0.03 

0.003 

Pressure Ratio 

90 

120 

600 

1200 

Temperature Ratio 

620/300 

720/300 

800/300 

900/300 


Using these values as input parameters the critical operating envelop 
can be established and the results are summarized in Figure 58 . Depending 
upon the initial operating conditions, a critical time between 0.1 ^sec 
and 10 msec is required. 
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In summary, the hot spot mechanism outlined above is capable 
of providing an explanation for all the major phenomena observed in 
the test programs. The pressure dependence (i.e., high altitude 
effect) . is attributed to the requirement of a high pressure ratio in 
order to exceed the minimum activation energy necessary for the de- 
composition of HN^ or other similar compounds. The characteristically 
high brisance of a relatively poor explosive is a result of the simul- 
taneous initiation of a number of localized "hot spots" contained 
throughout the residue. These "hot spots" maybe created either 
through a normal engine pressurization process during which the temper- 
atures reached are low but the time available for coupling the energy 
to the primary explosive is long, or alternatively through a gas phase 
detonation which results in a higher temperature for a much shorter 
time period. Additional experimental shock tube studies to further 
define which of these two possibilities is most important should be 
performed . 

6.4 Nonequilibrium Processes Associated with Rocket Ignition 
at Low Pressure 

6.4.1 Chemical Disequilibrium • 

The experimental results described in SectionsIV-V indicate 
that the assumption of equilibrium conditions is probably not valid 
for pulse-mode rocket engines operated at space ambient conditions. 

The calculated high temperature chemical composition for an 
O/F of 2.0 for an MMH/NTO propellant system show sizable concentrations 
of OH, H 2 O and HCO and only small quantities of NO 2 . Time integrated 
emission spectra, however, show strong band structures which are best 
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identified with a nitrogen dioxide NO 2 fluorescence, an HNO emission, 
and weak band of CH, EH, N 2 and other species. The absence of HCO or 
HCN bands in the emission is probably significant since for both species, 
the transition probabilities are high which should cause relatively 
strong spectral emission characteristics, even though their absolute 
concentration levels are small. The absence of these bands indicated 
that the reactions do not proceed to completion and that the system is 
not in chemical equilibrium. 

The conclusion that the system is not in chemical equilibrium is 
also supported by the discontinuous emission experimentally observed. 

For example, the emission record of sodium lines shows a series of 
discrete ignitions which are in fact closely spaced temperature pulses. 

6.4.2. Nonequilibrium Effects Caused by Disintegration of 
Propellant Jets at Low Pressures 

Ignition ’delay models describing the physical processes governing 
the kinetics of droplet evaporation and chemical processes governing 
the rates of chemical reaction prior to ignition have been applied to 
the problems of hypergolic ignition in space ambient engines. Mathe- 
matically in one model (Seamans), the propellant spray is represented 
by a three drop size array. The three drop sizes are obtained by 
applying a logarithmiconormal distribution to a given mass median 
drop size and then selecting three radii such that a prescribed weight 
flow distribution is satisfied. In this model both fuel and oxidizer 
drop sizes are assumed to follow the same logarithmiconormal distri- 
bution function. 
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The assumptions for this model were experimentally tested by 
observing the development and growth of propellant sprays under vacuum 
startup conditions. A photographic record of a typical high altitude 
bipropellant spray has already been given in Figure 36 . 

The oxidizer (^0^) fans out in a very fine mist, whereas the fuel 
(^H^) retains its jet integrity even at very low pressure. 

Based on this photographic evidence, several conclusions can 
be drawn. First, mixing in a one-on-one unlike injector used in this 
program is poor under vacuum startup conditions. In effect, the re- 
actants consist of a liquid column of hydrazine burning in a finely 
divided nitrogen tetroxide aerosol. Since a typical firing pulse of 
an attitude control engine is 20 msec in duration, a signifi- 

cant proportion of the fuel is able to strike a cold surface prior to 
ignition and this may be one of the reasons that large amounts of 
unreacted hydrazine were found in the residue mixtures. A second 
conclusion is that current theoretical models predicting the rate of 
initial pressurization of ^H^/NTO rocket engines must use different 
particle-size distribution functions for the fuel and oxidizer com- 
ponents if they are to be applied to problems involving vacuum startup. 

A series of calculations were performed in order to outline the 
effects of temperature j droplet size, and mass flux on rate of evapor- 
ation of rocket propellants. For these calculations it was assumed 

that the rate of instantaneous evaporation for a spherical liquid drop 

27 

can be calculated using the relation 
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(6. 4. 2-1) 


W a = 6 k n D 7 — 

a l-x a 

o 

where 

D is the diameter in cm 
x a# is the vapor pressure of A at the surface droplet 
x a* is the pressure of A at a long distance 


For forced convection k„ m is given by- 


Cf Dab 

[*2.0 + 0.6( DVooPf V 

2. vl/3 

( c -; ) 

D 

\ M* / 

V k / f J 


in which the subscript I denotes the properties of the fluid at 
the "film temperature" and "film composition" evaluated by 


T f ~ ^ ( T ° + T <» ) 

X af = ji(* a o ** a oo ) 


( 6 . 4 . 2 - 3 ) 


( 6 . 4 . 2 - 4 ) 


For high transfer rates, 0 , the correction coefficient is cal- 
culated from 


R • ( 6 . 4 . 2 - 5 ) 


where R denotes generalized flux ratio for momentum, energy and 
molar fluxes . 

Limiting the argument to high mass transfer rates becomes 


0 


j bn, ( Rqb+ 1 ) 

Rab 


(6. 4. 2-6) 
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in which R ab is given by 


Rab = 


Xqp ~ Xqoo 

N 


CO 


Nao + N b0 


- x, 


ao 


(6. 4. 2-7) 


where 

N „ is the molar flux of A at the surface 
N b0 i s 'the molar flux at B at the surface 
Under conditions where there is no mass transfer of B at the sur- 
face these equations reduce to 


=^tfRcb + Vnr 2 k xm 


(6. 4. 2-8) 


Physical parameters necessary for the evaluation of were 

either taken from literature values or were calculated for binary 
system. The value for Dqb ^was calculated using the expression: 


2.265 x 10 5 \/T( 1 /M a + 1 /M b ) 


'ab 


ab 


n 


( 6 . 4. 2 - 9 ) 


D ab 


where 

°oh is the Lennard Jones characteristic diameter, 

S2 d k is the Stockmayer potential taken as a 
temperature dependent function. 

Based on this theoretical formulation a series of calculations 
were performed. 

Typical results of this calculation are given for and in 

_2 

Figure 59 for a 10 mm drop size. These results provide an esti- 
mate of the initial mass flux from the spherical drop assuming the 
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VAPORIZATION moles /second 




ambient pressure given on the abscissa. The relative rates of 
evaporation differ by about two orders of magnitude at 100 psia 
and by about one order of magnitude at the lower start pressure. 

In addition to the pressure effect the effect of temperature is 
much stronger for hydrazine than for nitrogen tetroxide. 

In order to determine whether an improvement in the rate of 
dispersion of the hydrazine could be effected by increasing the 
temperature a series of experimental tests were performed. The 
temperature of the hydrazine was increased from 'room temperature 
to 120°F and l60°F and the hot hydrazine injected into the vacuum 
chamber . 

At 120°F the calculated instantaneous rate of evaporation of a 
droplet of hydrazine was estimated to be approximately the same as 
that of nitrogen tetroxide at room temperature and breakup of the 
fuel into a fuel-spray was followed using high speed photography. 

There was no observable effect caused by the increased temperature. 

The conclusion was reached that an increase in the injection tempera- 
ture of hydrazine to a temperature at which the initial rate of vapori- 
zation was calculated to be the same as the rate of vaporization for 
KTO under normal conditions did not result in an increased dispersion 
of the liquid stream. The theory used for this purpose is believed to 
represent a reasonably valid approximation and should be capable' of 
describing qualitative effects. The fact that this treatment, based 
on equilibrium vapor pressure and gas transport arguments apparently 
does not describe even qualitatively the experimental observations, 
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tends to indicate that other factors, e.g., dissolved gases or sur- 
face tension effects caused by trace impurities, may be responsible 
for the spray breakup. It is noted that relatively few studies have 
been performed outlining these effects at very low pressures. 

6.5 Flash Photolytic Effects 

An experimental problem in this program was to develop a 
suitable spectroscopic method which would allow the course of re- 
action to be followed by observing "simultaneously" time integrated 
emission-absorption spectra. An attempt was made to solve this prob- 
lem by the use of a short duration flash technique in which a very 
intense pulsed light source was used to backlight the rocket test 
engine at preselected frequency rates. A potential criticism of this 
approach was the possibility that the light pulse might flash photo- 
lyze and/or sensitize the reaction mixture bj' creating high free 
radical concentration. For any sensitization process two mechanisms 
are operative. One is an increase in the temperature of the reaction 

mixture caused by the thermalization of "hot" reaction products. The 

i 

second is an increase in the concentration of intermediate species 
with the resultant change in reaction velocity. In either case, if 
the system perturbation is sufficiently high, so that the system is 
brought to a point above the explosion limits, then the measurement 
methods would not be passive with respect to the reaction system. An 
estimate of the effects of photochemical activation was required, and 
an outline of the general theoretical approach necessary to evaluate 
these effects is presented in the following discussion. 
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The general’ approach used was to examine the time dependent 
behavior of a reaction system subjected to either short-duration or 
steady-state radiation fluxes. Because of the extreme complexity 
of this problem for the reaction mixtures present in a rocket chamber, 
preliminary calculations were restricted to each of the individual 
components. Fuel components which were considered include hydrazine 
(NgH^), monomethyl hydrazine (MMH), and unsymmetrical dimethyl- 
hydrazine (UDMH). Oxidizer components consisted of dinitrogen 
tetroxide ( ) and NO. 

A possible kinetic path for the thermal decomposition of 
hydrazine is outlined below: 


Initiation 


k, 

N 2 H 4 2NH 2 

(6.5-1) 

Propagation 

N 2 H 4 + NHg NH 3 + N 2 H 3 

(6.5-2) 

N^ N 2 + + H 

(6.5-3) 

N 2 H^ + H NH 3 + NH 2 

(6.5-4) 

Chain Branching 

N 2 H 3 NH + NHg 

(6.5-5) 

N 2 H^ + NH — N^ + N 2 H 3 

( 6 . 5 - 6 ) 

Termination 

k 7 

NH 2 + n 2 h 3 -4 nh 3 + n 2 + h 2 

(6.5-7) 

N 2 H 3 + NgHj 2NH 3 + N 2 

( 6 . 5 - 8 ) 

N 2 H 3 + H N 2 + 2H2 

(6.5-9) 

NHg + NHg N 2 H 4 

( 6 . 5 - 10 ) 
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Since this reaction sequence is complex, it is useful to adopt a 

28 

simplified sequence similar to that proposed by Berlad 
This includes an additional initiation step associated with a 
photochemical activation mechanism. 


Initiation 


N 2 H 4^ L *‘ 2NH 2 

(6.5-11) 

N 2 H 4 * 2m 2 

(6.5-12) 

Propagation 


; n 2 h 4 + KHg 

(6.5-13) 

Chain Branching 


N 2 H 3 — , HH ’ + NH 2 

(6.5-14) 

N 2 H4 + NH — *NH 2 + N 2 H 3 

(6.5-15) 

Termination 


N 2 H 3 + M2-»NH 3 + W 2 + Hg 

( 6 . 5 - 16 ) 

NH 2 + NH 2 -*-N 2 H4 

(6.5-17) 


A similar scheme was used for the decomposition of UDMH and MMH 
with appropriate modification of specific rate constants and 
energy factors . 

. The equations can be further simplified to a general form for 
the general series of hydrazine-based fuels. 


Photochemical 

Initiation 

F — ► A + NHg 

(6.5-18) 

Collisional Initiation 

F — *• A + KHg 

(6.5-19) 

Propagation 


F + A — ► P 1 + B 

(6.5-20) 

Propagation 

F + 

NHg — + B 

(6.5-21) 

Termination 

B + 

NH 2 — > Pp + n 2 + h 2 

(6.5-22) 
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In this simplified format the tabulation of species F, A, B and 
P 1 for these fuels is given by: 


Species 

Hydrazine 

MMH 

UDMH 

F 

n 2 H4 

CHjHNHHg 

(CH 3 ) 2 MH2 

A 

NHg 

CH^HN 

(CH^gN 

B 

"2=3 

CH^HN - NH 2 

CH^CHg - NHg 

P 1 

“3 

CH^HgN 

(CH^gHN 


For hydrazine, the sequence of reactions can be written: 

N 2 H 4”*' 2NH 2 (6.5-23) 

+ NH 2 — ► NH^ + NgHg (6.5-24) 

+ NH — ►NHg + (6.5-25) 

NgH- + NHg — ►HH 3 + N 2 + Hg (6.5-26) 

Using a steady state treatment for the intermediate radical 
species NH 2 , N 2 H 3 and NH, we solve to find the steady state con- 
centration of NH2 formed by thermal activation processes: 

[ ira 2] ss = T7 + T7 (6.5-27) 

If we include a photochemical activation mechanism given by 

f N 2 H 4 — + 2m 2 (6.5-28) 

then the photochemical steady state concentration for NHg is 
given by 

[my -- tr*-E7 (6 - 5 ‘ 29) 

A combined thermal and photochemical activation mechanism then .. 


k v k, k 5 

rr + + 

k 2 k- k 7 


1 46 


yields for NHg 
[NHg] = 


( 6 . 5 - 30 ) 



A further modification for collisional activation was 


developed to yield 

W + 4H m] (6 - 5 ' 31) 

where M represents the activating third body. The explicit use 
of this latter formulation aids in the formulation of pressure 
dependencies for the reaction sequence. Since the vacuum ignition 
involves a reaction which occurs throughout a wide pressure regime 
the rate constants have to be examined to determine whether they 
are themselves functions of pressure over the range of the ignition 


process. 

The effect of flash photolysis was determined by comparing 
the rate of decomposition for individual components, which is 
associated with the energy density flux of the light source, to 
that which is expected for the normal thermal decomposition. 

This was evaluated by examination of the extinction coefficients 
for specific components and intermediate species for selected 
flash intensities. 

'Calculations were performed to assess the effects of possible ■ 
photochemical activation on a change in reaction mechanism. Results 
of this calculation are similar to those of Berlad, corresponding to 
a gas density for at 60 psia and 300°K. Four plots are presented: 

(l) the reaction rate vs. time (Figure 60 ); (2) relative hydrazine 

concentration vs . time (Figure 6l); (3) HH vs. time (Figure 62); 

and (4) RHg.vs. time (Figure 63). As can be determined from the 
plot of reaction rate vs, time, an induction time is observed 
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corresponding to a time delay of about 150 microseconds. During this 
interval, the hydrazine concentration is reduced by about 15 percent 
and a corresponding buildup of reactive intermediates occurs. 

The major conclusion of this analysis was that there should 
be no appreciable perturbation of the ignition process caused by the 
absorption of radiant energy by hydrazine. However, it should be 
noted that it may be possible to alter the reaction mechanism by 
using somewhat higher energy levels over a broader geometric cross 
section. Such a technique could provide a way of significantly 
reducing ignition delay times and thereby reduce the amount of resi- 
dues formed during pulse-mode operation. 
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6.5*1 Comparison of Gas Parameters and Ignition Time 
for N_H, and MMH 

2 4 

2 

Tbe Seamans program was used to calculate gas temperature, and fuel 
and oxidizer vapor pressure in the reaction chamber. Results for N^H^/EITO 
and MMH/NTO are presented in Figure 64. The calculations indicate that 
the ignition time for is about 35% of that for MMH. 

Variation of several program parameters produced no appreciable change 
in the program results. The parameters varied were: fuel droplet temper- 

ature, fuel evaporation coefficient, activation energy for ignition (this 
affects ignition delay time only), activation energy for intermediate, 
and maximum and minimum values of drop wall heat transfer factor. 
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Figure 64. Calculated Gas Temperature and Fuel and Oxidizer Vapor 
Pressure for N^/NTO and MMH/NTO 


152 


GAS TEMPERATURE - °R 


6.6 Relative Activation Energies for NpHjj. , MMH and UDMH 

The "hard start" pressure transients associated with the ignition of 
hydrazine based fuels decrease in the order > MMH> UDMH . This can be 
correlated to the decrease in rate of decomposition found by Eberstein and 
Glassman"^ in temperature ranges between 800 and 1220°K where the first 
order rate of reaction took the ordering UDMH >MMH> N^H^. Extrapolation of 
their data shows that there is a crossover between the rate of reaction of 
^Hj^ and MMH at about 750°K and a crossover for MMH and UDMH at about 1100°K. 
Presumably the hard start problem is related in some way to these rates of 
decomposition and where the rate of decomposition is fast the ignition 
proceeds more smoothly. Experimental results of Eberstein and Glassman 
show that the activation energy required for decomposition had the order 
MMH> NpH^UDMH. Assuming that the rate limiting step corresponds to the 
initial N-N bond dissociation then the predicted ordering of the relative 
activation energies would be UDMH> MMH>N 2 H^. 

29 

Electron impact studies by Dibeler et al , show that the N-N bond 
is strengthened by the addition of methyl radicals and that the (N-N) bond 
strength of N^H^, MMH, and UDMH are 60, 67, and 72 Kcals/mole respectively. 
iyro lytic of hydrazine using a toluene carrier (Kerr et al^®) showed the 
opposite ordering for the N-N bond of 57.1, 51-9 and 49.6 Kcals/mole 
respectively. Kerr's data is in agreement with the kinetic results of 
Eberstein and Glassman and with the observation in the rocket propulsion 
field. Dibeler' s data would predict the opposite trend unless a chain 
branching mechanism involving a methyl radical would cause a change in 
kinetic mechanism through this homologous series. Ebersten and Glassman 
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and most of the subsequent workers in the propulsion field have developed 
a model for the decomposition of hydrazine fuels in attempts to explain 
the rough start in terms of a branching model in part because it is then 
possible to reconcile the results of Kerr and Dibeler. 

However, it is not necessary to resort to a branching chain mechanism. 
The problem arises by the use of the simple Arrhenius model describing 
reaction rates in which 

k = A exp (-E/RT) 

In this relation the reaction velocity, k, is proportional to an Arrhenius 
preexpotential factor A and a quantity exp (-E/RT) where E is the activa- 
tion energy and R and T are the gas constant and absolute temperature 
respectively. Using the Arrhenius form of kinetic behavior and correlating 
experimental results with theory experimental results are interpreted in 
terms of an energy of activation and a preexpotential factor, A. Both of 
these are usually assumed constant for engineering calculations. Using 
more detailed theories of kinetics, as derived by methods of either clas- 
sical or quantum statistical mechanics , it can be shown that the experimental 
activation energy E is related to the probability that the molecule contains 
a total energy E' distributed among y-square terms. Since each effective 
mode of vibration contributed two square terms the maximum number possible 
for a nonlinear x-atom molecule y has the value 

3n - 6/2 

In practice, the number of effective square terms are chosen relative to 
experimental data and y becomes a somewhat arbitrary parameter related to 
the complexity of the molecule. 'The implication of this is that the 
experimentally observed activation energy is not only a function of bond 
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strength hut also is dependent on the ability of the molecule to exchange 
stored energy into the bond. This means that all other things being equal 
a more complicated molecule will exhibit a lower activation energy and 
that the experimental rates should follow the order UDMH> MMH> NgH^ which 
is, in fact, observed. This type of argument can be us.ed and kinetic 
models developed which do not necessarily have to resort to chain branching 
mechanisms . 

The reaction velocities for decomposition also exhibit a crossover 
at about 750°K for MMH and and at about 1250°K for MMH and UDMH. A 

possible explanation for this crossover is that the number of effective 
oscillators which exchange energy with the breaking bond increase with 
increasing temperature and that the number of possible interacting cross 
terms increase with increasing molecular complexity. The detailed con- 
tributions of interacting nonorthogonal cross terms are not predicted by 
any present theory. Therefore, it would be expected that the rate of de- 
composition of more complex molecules would increase more rapidly than less 
complex molecules. Thus, it would be expected that MMH reacts less rapidly 
at lower temperatures than and that a crossover at 750°K would occur. 

In such a detailed analysis it would be expected that at low temperatures 
corresponding to the fuel rich region surrounding a fuel droplet the rate 
of reaction of N^H^ would exceed that of UDMH and produce substantial 
quantities of an unstable chemical intermediate. These effects’ would be 
greatly magnified for pulse mode operations. 

One result of the above argument is that it demonstrates how difficult 
it is to develop a kinetic model for the decomposition of hydrazine propel- 
lants without a much greater understanding of the fundamental chemical 
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reaction rates. The use of the reaction mechanism which may be valid at 
1000°K may not necessarily describe the expected rates of reaction at low 
temperatures surrounding the fuel droplet, and the use of the kinetic 
mechanisms describing the combustion processes at higher temperatures may 
or may not be valid. The use of global kinetics where a constant activation 
energy is assumed or experimentally evaluated may not allow extrapolation to 
more or less complicated molecules. Since with very few exceptions most 
engineering approaches use either a first or second order kinetics scheme 
to simplify mathematical solution, the results from these engineering calcu- 
lations are a priori extremely questionable. At the present time, there is 
no conclusive evidence that the kinetic mechanism usually assumed for the 
decomposition of hydrazine fuels correctly describes the process since it 
is possible to explain the observed chemical stoichiometry and reaction 
orders using an alternative mechanism involving the unimolecular decomposi- 
tion of N^H^. Such a decomposition would result in the production of meta- 
stable azide derivatives in both the gas and liquid phases . 
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VII. CONCLUSIONS 


Based on the results of this program, several conclusions 
can he formulated. 

1. The catastrophic overpressures observed in space simu- 
lated, pulse-mode, rocket engines using nitrogen tetroxide/hydrazine 
propellants are a result of a liquid phase detonation in the un- 
burned propellant residues formed during previous firing. The 
mechanism of detonation involves one or more "hot spots” generated 
by an adiabatic compression of gas bubbles contained within the 
residues. These "hot spots" undergo an exothermic chemical reaction 
forming a localized high energy region in the propellant residues 
and then transfer sufficient heat to cause explosive initiation in 
the residues. 

2. A chemical characterization of the propellant residues 
indicates that they consist principally of a mixture of fuel and 
fuel nitrates. In addition, a number of other chemical species 
are also present. In the liquid phase these include hydrogen 
azide, triazene, hyponitrous acid and tetrazene. In the gas phase, 
hydrogen azide, triazene together with other unidentified species 
are present. 

3. The usual assumptions of chemical and physical equilibrium 
used in engineering design calculations do not apply during the ig- 
nition of space ambient rocket engines using NTO-Hydrazine propellants. 
Physical disequilibrium is caused by a flash vaporization of the 


157 



nitrogen tetroxide resulting in the formation of a finely divided aerosol 
whereas, the fuel stream remains as an integral highly turbulent jet. 

These effects result in strong chemical disequilibrium which allows the 
formation of sensitive explosive intermediates probably via a heterogeneous 
reaction between NO^ and liquid fuel. 

4 . The use of time resolved spectrophometric techniques can provide 
information concerning intermediate chemical species after ignition has 
occurred. 

5. The use of time resolved spectrochemical techniques could not 
provide information concerning intermediate chemical species in the pre- 
ignition time regime. 

In summary, the major work effort of this research program has been 
devoted to finding a general mechanism which can cause the initiation of a 
propellant residue. The mechanism suggested is capable of explaining all 
the phenomena thus far observed. In particular it provides a unique 
explanation for observed high brisance of the explosive material and pre- 
dicts the observed high altitude dependency. 
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APPENDIX A 

Specific Reaction Rate For Nitrogen Tetroxide, Hydrazine, and Their 

Decomposition Products 

Figures 1 through 25 show the variation of specific reaction rate 

constant, k, with temperature for selected gas phase reactions involving 

nitrogen tetroxide and hydrazine. The units of k are: sec ^ for uni- 

3 -1 -1 

molecular reaction, cm - mole - sec for bimolecular reaction, and 

6 - 2-1 * 
cm - mole - sec for termolecular reaction. 

In many cases, experimental variations of several orders of magnitude 

sure reported. Since it was beyond the scope of this program to critically 

evaluate each of the reported references, the variations are reported 

without comment. It is perhaps instructive to note the wide variations 

as an example of the lack of precision and accuracy in the field of gas 

phase kinetics studies. 

A list of references from which the reaction rate data were obtained 
is given on pages A27 - A33« This is followed by a table listing those 
references that were used for the reactions in each figure. 
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Calibration Spectra In Visible Region 

-12. Argon Discharge Tube Spectra 
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Test #14. . Monopiethyl Hydrazine/Nitrogen Tetroxide. 
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Test #15. 


Monomethyl -Hydrazine/mtrogen Tetroxide - Short Oxider Lead.. 

a. Pressure and visible emission spectra vs. time. Time 
base 500 microsecond per cm. Long wave length visible 
emission spectra. • 

"k* Pressure and visible emission spectra vs. time. Time 
base 500 microsecond per cm. Short wave length visible 
emission spectra. 

c. Visible emission spectra vs. time. Time base 500 micro- 
second per cm. Upper trace short wave length visible 
emission. Lower trace long wave length visible emission. 

d. Visible emission spectra vs. time. Time base 100 micro- 

seconds per cm. Upper trace short wave length visible 
emission spectra. Lower trace long wave length visible 
emission spectra. ~ 

e. Pressure and long, wave length visible emission spectra 
vs. time. 
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Test # 16 . Monomethyl Hydrazine/Nitrogen Tetroxide - Fuel Lead. 

a. Pressure and short wave length visible emission spectra 

vs. time. Time base 500 microseconds per cm. 

b. Visible emission spectra vs. time. Time base 500 micro- 
seconds per cm. Upper trace short wave length. Lower 
trace long wave length. 

c. Pressure and visible emission spectra vs. time. Time 
base 500 microseconds per cm. Upper trace pressure. 
Lower trace long wave length visible emission spectra. 

d. Visible emission spectra vs. time. Time base 100 micro- 
seconds per cm. Upper trace short wave length emission 
spectra. Lower trace long wave length emission spectra. 

e. Pressure and visible emission spectra vs', time. 
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Test #18. 


* Monomethyl Hydrazine/witrogen Tetroxide. Oxider Lead. 

a. Chamber pressure and visible short wave length emission 
spectra vs. time - time base 100 microsecond per cm. 

b. Chamber pressure and visible emission spectra 

. > 

(4000 - 6900$) vs. time - time base 100 microsecond 
per cm. unless otherwise indicated. 

c. Visible emission spectra vs. time - time base 500 
microseconds per cm. Upper trace.short wave length. 
Lower- trace long wave length. 

d. Visible emission spectra vs. time - time base 100 
microseconds per cm. - Upper trace short wave length. 
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Lower trace long wave length. 
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Test #23. Unsymmetrical Dimethylh'ydrazine / Nitrogen Tetroxide 
Fuel Lead. 

a. Short and long wave length Emission Spectra in 

' visible region - time base 500 microseconds per 

. > 

cm. 

b. Short and long wave length Emission Spectra in 
visible region - time base 100 microsecond per 
cm - 250 microsecond dead time between oscillo- 
graph tractings. 
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Test #30. Hydrazine/Nitrogen Tetroxide, Short. Oxider Lead. 

a. ‘Chamber pressure and propellant flows vs. time. 

Time base 20 milliseconds per cm. 

b. Chamber pressure and IR emission spectra vs. time. 
Time base 500 microseconds per cm. 

c. Chamber pressure and IR emission spectra vs. time. 
Time base 100 microseconds per cm. 250 microseconds 
dead time between oscilloscope records. 
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APPENDIX D 


Spectral Characteristics of Selected 
Major Combustion Products 
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APPENDIX D. SPECTRAL CHARACTERISTICS OF SELECTED 
MAJOR COMBUSTION PRODUCTS 
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APPENDIX D (Continued) 


Swan Mulliken 
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4075 

6784 





4267 . 

7231 





5145 

7236 





43130. 





2816 

Red 

Main 


Tail 



63550 

3584 

3883 

3405 

3629 

2164 

5502Q 

3586 


3433 

3880 

2069 

WTr 

3590 

4i8i 

3465 

3910 

I8l4 

8067 

3^55 

4197 

3603 

39^5 


9140R 

3862 



3985 


9393 

3871 






l64i 

1608 

1107 



APPENDIX D . (Continued) 

Species Visible (%) IR (cm" 1 ) 

CN" 2250-2050 


C 2 N 2 

U. V. 


2322 





2149 


0 

0 

+ 


•• . . ' - 

2214 





1734 




■ . • - . 

- 1562 


H 3 BO 3 



3250 

1060 




3150 

881 




l44o 

652 



1 ■ 

. 1185 

544 

HCR 

7880a 


3311 



7912a 


2097 





712 


HUGO 


. 

3531 

; 797 




2274 

670 




1327" 

572 

HRO 



3450 





trans 

cis 

HNO 

3417a 


3426 

3590 

u 

3539a 


i64o 

1696 


3545a 


1292 

1260 


3680a 


856 

794 




637 

598 





543 

h 2°2 

3700a 


3610 

1262 




3417 

890 




2870 



Flames 

Chemiluminescence 



NH 

3360Q 

3035R 

3300 



3370Q 

3042Q 




324or 

32530 

3610R 

36270 



APPENDIX D (Continued) 


Species 


Visible (a) 


IE (cm -1 ) 

nh 2 

4723 

5707 





5166 

5708 





5525 

• 6333 

. 




5705 

’ ■* 




NH 




N l4 

N 15 

J 




3448 

3335 



- 


34l4 

1625 





3337 

961 



* 


1627 

926 





968 





• 

950 




* 


932 


n 3 h 




3336 

1150 




2140 

672 





1274 

522 

NO 

3886 

* 


2374 

1712 


3801 



2371 

1380 





2347 

1300 





2327 

1262 





1904 14 

I876 - N ;■ c 

1843 -n 5 

1170 





1155 



- 


1038 

■ - 




1748 

1037. 

1020 


stretch: 2400-900 
monomer : 1883 
cis-dimer: 1862, 
1768 

trans -dimer: 1740 

NO 6000 2220 


d4 


-APPENDIX D (Continued) 


Species 


no 2 


no 2 + 


n 2 o 


Visible (8) 


4350a 

4390a 

mk 

Woa, 

4630a 


IR (cm" 1 ) 




1618 

485 

1580 

1318 

482 

1306 

808 

429 

740 

750 

381 


672: 

370 



0“N-0 stretch: 900-400 


2360 

1400 

667 


Continuous absorption through 
U.V. to 30658 


/W 


n 15 h 15 o 


2224 

2178 

1286 

1281 

. 589 

576 

f l 5 n 1 ^o 

n 15 n 15 o 

2203 

2156 

1271 

1266 

586 

572 


D5 


APPENDIX E 


Visible Coincidence Table - Minor Species 










'kp 


f^O 


8 


APPENDIX E VISIBLE COINCIDENCE TABLE 
MINOR SPECIES 


3035R 

- 

NH 


3910 

- 

CN 


5586 

- C 2 

30420 

• 

NH 




+ 


5636 


3065 a 

- . 

N„0 


3921 

- 

C 


- o 2 



CL 


3945 

- 

CN 




324or 

• 

NH 


3985 

- 

CN 


5705 

- NHp 

32530 

• 

NH 




_1_ 




3273 

- 

B ? 


4075 

- 

C 


5707 

- nh 2 





4l68 

- 

CN 

- 



3360Q 

* 

NH 


4l8l 

- 

CN 


5708 

- nh 2 

3370Q 


NH 


4197 

- 

CN 



3386 

- 

NO 


4216 

- 

CN 


5794 

- c 

3^05 

- 

CN 







I 

3417a 

- 

HNO 


4267 

- 

C + 


6000 

- NO 



2 


43130, 

- 

CH 


6333 

- NHp 

3433 

- 

CN 


4350 

- 

NOo 



c. 

3465 

- 

CN 




u 


63550 

- CN 

3539a 

• 

hno 2 


4390a 

- 

N0 2 


65020 - 

•- CN 

3545a 

• 

hno 2 


4448a 

-• 

CM 

0 


6578 

- C + 

3584 

_ 

NO,CN 


4480a 

— 

NOp 


8583 

- c + 

3586 

- 

CN . 








3590 

- 

CN 


4630a 

- 

no 2 


6784 

- c + 

3603 

- 

CN 






7119 

- c 

3610R 

- 

NH 


• 4698 

- 

Cg 



4 . 

36270 

- 

NH 






7231 

- c 

3629 

- 

CN 


4715 

- 

C ? 



+ 

3680 a 

- 

HNO„ 




dL 


7236 

- c 

- 


2 


4723 

- 

m 0 


7874 

- CN 

3700a 

~ 

HpOp 


4737 ’ 


c. 


7880 a 

- HCN 





- 

c ? 


7912a 

- HCN 

3768R 

- 

BH + 




C. 


8067 

- CN 



+ 

- 

5129 

- 

C 2 


914or 

- CN 

3792Q 


BH 




+ 


9393 

- CN 

3801 


NO 


5145 

- 

c 




"I • ' i ^ ' 

- 







(a) - 

absorptii 

3803R 

- 

BH + 


5165 

- 

Cp 


(0) - 

Q branch 

3855 

- 

CN 




c. 


(R) - 

R branch 

3862 

- 

CN 


5166 

- 

NHp . 



3871 

~ 

CN 




CL 




3876 

- 

C+ 


5525 

- 

nh 2 




3880 

- 

CN 








3883 

- 

CN 


5541 

-• 

c 2 





El 



3610 - h 2 o 2 

3590 - EN0 2 

3531 -HNCO 
3450 - HUO 
3448 - m 3 

3426 - M0 2 

3417 - h 2 o 2 

34i4 - uh 3 

3337 - i® 3 

3336 - N^ 

3335 - nh 3 

• 3311 - HCN 
3300 - NH 
3250 - H 3 B0 3 

3150 - H 3 B0 3 

3000 - b(oh) 2 

2976 - BH 3 

2870 - h 2 0 2 

2816 - CH 
2650 - EH 2 

2500 b(oh) 2 

2430 - bh 2 

2384 - bh 3 

2374 - NO 
2371 - NO 

2360 - N0 2 + 

2347 - NO 
2327 - W0 


APPENDIX E (Continued) 


INFRARED COINCIDENCE TABLE 

2322 - c 2 n 2 

2274 - HNCO 
2270 - B 10 H 

,, 4 

2264 - B 11 ^ 

2250 - B 10 ^ 

2244 - B 1 ^ 

?224 - N 2 0 

2214 - C0 + 

2203 - N O 

1 ([ 

2178 - n 2 o , 

2164 - CN 

2156 - n 2 o 
2149 - C 2 N 2 
2l40 - n h 
•2097 - HCN 

2069 - CN 
1904 - no 
1876 - no 
1856 - c 2 

1843 - NO 

1830 -c 2 

18 14 - CN 
1809 -C 2 

1788 - c 2 

1765 - BHj 

1748 - NO 

1734 - C0 + 

1712 - NO 

I696 - MQ, 


MINOR SPECIES 

i64o - hno 2 

.1627 - nh 3 

1625 - NH 3 

1618 - N0 2 

1608 - c 2 

1580 - N0 2 

1562 - co + 
i44o - h 3 bo '■ 

1400 - NO + 

2 

1380 - NO . 

1327- HNCO 
1322 - B0 2 

1318 - no 2 
1306 - no 2 

-3 

1300 - NO, B0 3 ,B(0H) 2 

1292 - MO 

2 

1286 - n 2 o 


1281 

1274 

1271 

1262 

1260 

1210 

1208 

1185 

1170 

1155 

1150 

1107 

1100 

' 1093 

1080 

1070 

106l 

1060 

1038 

1037 

1020 

1000 

.968 

961 

950 

9^7 


APPENDIX E. (Continued) 


n 2 0 

932 

- nh 3 

586 

- 

n 2 o 

V 

926 

-nh 3 

576 

- 

n 2 o • 

n 2 0 

890 

- ^ 

576 

- 

N 2 0 . 

h 2 o 2 , no 

- 881 

- H 3 B0 3 

572 

- 

HNCO, ] 

hno 2 , n 2 o 

856 

- HN0 2 

544 

- 

H 3 B0 3 

B n H 4 

840 

- BH 3 

543 

- 

hno 2 

B 10 H U 

808 

- N0 2 

533 

- 

b(oh ) 4 

H^BO^ 

802 

- bh 3 

522 

- 

N3H 

NO 

797 

- HNCO 

485 

- 

NOo 

NO 

794 

- hno 9 



d 

B(OH), 

754 

d 

- B(OH), 

482 

— 

N0 2 

C 2 



464 

- 

BOo 

d 

750 

- NOo B0o" 3 ,B(0H)o“ 



d 

b(oh ) 2 



460 

- 

B(0H)o 

b 1 \ 

749 

* B(0H) 4 

450 


d 

B(0H)o 

B 1 \ 

747 

- B(0H) 4 " 

429 

_ 

d 

N0 2 


740 

- N0 o 




B0 2 


d 

381 

- 

no 2 


712 

- HCN . 




B 2 

672 

- N 3 H,N0 2 

379 

- 

b(0H ) 4 

H3BO3 

670 

- HNCO 

370 

- 

no 2 

NO 

667 

- no 2 + 




NO 






NO 

652 

- H a B0 Q 




B(OH) 


3 3 




d 

650 

- BO q " 3 




NH 0 






3 

637 

- HNOo 




NH, 


C— 




3 

598 

-■HN0 o 




B(OH) rNH 0 

589 

d 





- N 5 0 




b(oh ) 4 


d 





E3 


